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**  is.  (continued)  probability  of  crack  exceedance,  cumulative  distribution 
of  service  time. 


19.  (continued)  methods  for  determining  the  initial  fatigue  quality  of 
fastener  holes,  the  probability  of  distribution  of  service  time  to  reach  any 
specified  crack  size;  ^4)  provide  guidelines  and  design  data  for  implement¬ 
ing  the  durability  analysis  method  and  for  assisting  contractor  and  USAF 
personnel  in  complying  with  the  intent  of  the  durability  specifications  for 
metallic  airframes,  TA-c,  mc^ _ _ 

The  durability  analysis  methodology'^developed  provides  a  ^idurabil-ity?1  design 
tool^^-for  quantitatively  reflecting  durability  requirements  in  the  design 
process  and  for  making  design  tradeoffs.  The  methodology? accounts  for  the 
initial  fatigue  quality  variation,  crack  growth  damage  accumulation  in  a 
papulation  of  structural  details  (e.g.,  fastener  holes,  lugs,  fillets,  cut¬ 
outs,  etc.),  load  spectra  and  structural  properties.  . 

During  manufacturing  and  assembly,  flaws  of  various  types n(a.g. ,  scratches, 
burrs,  microscopic  imperfections,  etc.)  care  produced  in- structural  details,  . 
such  as  fastener  holes,  cutouts,  etc.  ^ThVTnTtXaI~f a tigue  qua  1  i ty  o faucir* 
details  is  represented  by  an  equivalent  initial  flaw  size  distribution 
(EIFSD) .  An  equivalent  initial  flaw  size  (EIFS)  is  an  artificial  crack  size 
which  results  in  an  actual  crack  size  at  an  actual  point  in  time  when  the 
initial  flaw  is  grown  forward.  ElFSs  are  determined  by  back-extrapolating 
fractographic  results.  C ede.')  J- 

The  EIFSD  is  the  "cornerstone"  fcr  the  durability  analysis.  Once  the  EIFSD 
has  been  determined,  a  two-segment  deterministic-stochastic  crack  growth 
approach  is  used  to  grow  the  EIFSD  forward  to  determine  either  the  probab¬ 
ility  of  crack  exceedance  at  any  service  time  and/or  the  cumulative  distri¬ 
bution  of  service  time  to  reach  any  specified  crack  size  x. .  The  predicted 
prob  ibility  of  crack  exceedance  can  be  used  to  estimate  statistically  the 
"extent  of  damage"  for  the  durability-critical  component.  These  include  the 
average  and  upper  bound  limit  of  the  "extent  of  damage"  for  selected  exceed¬ 
ance  probabilities. 

A  durability  analysis  methodology  has  been  developed.  Theoretically,  it  ap¬ 
plies  to  any  type  of  structural  detail  in  a  metallic  structure.  The  method¬ 
ology  has  been  demonstrated  for  straight-bore  and  countersunk  fastener  holes 
with  clearance-fit  fasteners  in  7475-T7351  aluminum.  A  comprehensive  demon¬ 
stration  has  been  conducted  for  coupon  specimens  and  for  a  full-scale  lower 
..  wing  skin  for  a  fighter  aircraft.  Applications  for  other  structural  details 
•  (e.g.,  cutouts,  fillets,  etc.)  need  to  be  investigated.  Also  fastener  holes 
:  with  interference-fit  fasteners,  cold-working,  etc.,  need  to  be  investigat- 
t  ed.  , 

The  data  for  defining  the  initial  fatigue  quality  for  different  materials 
and  structural  details  can  be  acquired  economically  and  timely  as  a  part  of 
the  Aircraft  Structural  Integrity  Program  (ASIP) .  Structural  details  in  se¬ 
lected  test  specimens  should  not  be  preflawed  so  that  baseline  data  can  be 
obtained  to  satisfy  the  data  requirements  for  initial  fatigue  quality,  dur- 
ability  and  damage  tolerance.  ■ 


_ Qnclasaifled _ __ 

3ICURITY  CLASSIFICATION  of  this  f*oj 


FOREWORD 


The  aacond  adition  of  thia  handbook  was  praparad  by  Gen- 
aral  Dynamics  (Fort  Worth  Diviaion)  and  by  United  Analysis 
Incorporated  (Vienna,  VA)  under  Phase  III  Taek  VIII  of  the 
"Advanced  Durability  Analysis"  program  (Air  Force  Contract 
F33615-84-C-3208)  for  the  Air  Force  Wright  Aeronautical  Lab¬ 
oratories  (AFWAL/FIBEC) .  Marge  E.  Artley  was  the  Air  Force 
Project  Engineer.  James  L.  Rudd  and  Dr.  Jack  W.  Lincoln 
(ASD/ENFS)  were  technical  advisors  for  the  program.  Dr. 
Sherrell  D.  Manning  of  the  General  Dynamics'  structural  Tech¬ 
nology  Staff  was  the  Program  Manager  and  co- investigator 
along  with  Dr.  Jann  N.  Yang  of  United  Analysis  Incorporated. 

The  following  reports  (AFWAL-TR-B6-3017)  were  also  pre¬ 
pared  under  the  "Advanced  Durability  Analysis"  program 


o  Volume  I  -  Analytical  Methods 


o  Volume  II  -  Analytical  Predictions,  Test  Results 

and  Analytical  Cci-i^lations 

o  Volume  III  -  Fractographic  Test  Data 

o  Volume  IV  -  Executive  Summary 


o  Volume  V  -  Durability  Analysis  Software  User's  Guide 


This  handbook  covers  work  accomplished  during  the  period 
October  1984  through  September  1987.  This  report  was  re¬ 
leased  for  publication  in  August  1988. 
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SECTION  I 

INTRODUCTION 

The  second  edition  of  the  Durability  Design  Handbook  for 
metallic  airframes  has  been  updated  to  also  cover  functional 
impairment  due  to  fuel  leaks  and  ligament  breakage.  In  the 
first  edition  [1]  the  durability  analysis  methods  were  limit¬ 
ed  to  functional  impairment  due  to  excessive  cracking  in  the 
small  crack  size  region  (e.g.,  cracks  <  0.10").  The  initial 
durability  analysis  method  for  the  small  crack  size  region 
has  been  extended  to  also  cover  large  through-the-thickness 
cracks  (e.g.,  0.50”  -  0.75”)  associated  with  fuel  leaks  and 
ligament  breakage. 

Objectives  of  the  handbook  are:  (l)  summarize  and  in¬ 
terpret  the  essential  U.  S.  Air  Force  Durability  Design  re¬ 
quirements  for  metallic  airframes,  (2)  provide  durability 
analysis  criteria  for  economic  life  and  durability  critical 
parts,  (3)  provide  state-of-the-art  durability  analysis  con¬ 
cepts  and  methods  for  determining  the  initial  fatigue  quality 
of  clearance-fit  fastener  holes,  the  probability  of  crack  ex¬ 
ceedance  at  any  service  time  and  the  cumulative  distribution 
of  service  time  to  reach  any  specified  crack  size,  (4)  provide 
guidelines  and  design  data  for  implementing  the  durability 
analysis  method  and  for  assisting  contractor  and  USAF  person¬ 
nel  in  complying  with  the  intent  of  the  durability  specifica¬ 
tions  for  metallic  airframes. 

The  methodology  accounts  for  the  initial  fatigue  quality 
variation  of  structural  details,  the  crack  growth  accumula¬ 
tion  for  a  population  of  structural  details  under  specified 
design  conditions  and  structural  properties.  Step-by-step 
procedures  are  provided. 

The  durability  analysis  approach  [2-6]  conceptually  des¬ 
cribed  in  Fig.  1-1,  reflects  a  probabilistic  approach,  a 
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aa:  Ref.  crack  size 

EIFSD  :  Equivalent  Initial  Flaw  Size  Distribution 


Figura  1-1.  Durability  Analysis  Approach. 


fracture  mechanics  philosophy  and  both  deterministic  and  sto¬ 
chastic  crack  growth  methods.  It  can  be  used  to  predict  the 
probability  of  crack  exceedance  at  any  service  time  and/or 
the  cumulative  distribution  of  service  time  to  reach  any 
given  crack  size.  The  methodology  applies  to  the  small  crack 
size  range  associated  with  excessive  cracking  (e.g.,  <  0.05") 
and  to  large  through-the-thickness  cracks  (e.g.,  0.50"  - 
0.75")  associated  with  fuel  leakage/ ligament  breakage. 

No  matter  what  form,  location  or  combination  the  as- 
manufactured  flaws  may  have  in  structural  details  (e.g., 
scratches,  burrs,  microscopic  imperfections,  etc.)  or  what¬ 
ever  the  source  of  fatigue  cracking  may  be,  a  practical  me¬ 
thod  of  representing  the  as-manufactured  condition  is  needed 
for  durability  analysis.  This  is  taken  care  of  by  the  equiv¬ 
alent  initial  flaw  size  (EIFS)  concept. 

An  equivalent  initial  flaw  is  an  artificial  crack  size 
which  results  in  an  actual  crack  size  at  an  actual  point  in 
time  when  the  initial  flaw  is  grown  forward.  It  is  determin¬ 
ed  by  back-extrapolating  fractographic  results.  It  has  the 
following  characteristics:  (1)  An  EIFS  is  an  artificial 
crack  assumed  to  represent  the  initial  fatigue  quality  of  a 
structural  detail  in  the  as -manufactured  condition  whatever 
the  source  of  fatigue  cracking  may  be,  (2)  it  has  no  direct 
relationship  to  actual  initial  flaws  in  fastener  holes  such 
as  scratches,  burrs,  microdefects,  etc.,  and  it  cannot  be 
verified  by  NDI,  (3)  it  has  a  universal  crack  shape  in  which 
the  crack  size  is  measured  in  the  direction  of  crack  propaga¬ 
tion,  (4)  EIFSs  are  in  a  fracture  mechanics  format  but  they 
are  not  subject  to  such  laws  or  limitations  as  the  "short 
crack  effect,"  (5)  it  depends  on  the  fractographic  data  used, 
the  fractographic  crack  size  range  used  for  the  back-extrapo¬ 
lation  and  the  crack  growth  rate  model  used,  (6)  it  must  be 
grown  forward  in  a  manner  consistent  with  the  basis  for  the 
EIFS,  and  (7)  EIFSs  are  not  unique  —  a  different  set  ^s 
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obtained  for  each  crack  growth  law  used  for  the  back-extrapo- 
1  lation. 

The  initial  fatigue  quality  (IFQ)  of  a  structural  detail 
(e.g.,  fastener  holes,  cutouts,  fillets,  lugs,  etc.)  is  re¬ 
presented  by  an  equivalent  initial  flaw  size  distribution 
(EIFSD) .  EIFS  is  treated  as  a  random  variable.  A  different 

E1FSD  should  be  established  for  different  types  of  structural 
details  (e.g. ,  see  Fig.  1-2)  reflected  in  the  durability  analy¬ 
sis  . 

The  durability  analysis  can  be  used  to  estimate  statis¬ 
tically  the  "extent  of  damage"  (such  as  mean  and  extreme  val¬ 
ues  in  a  durability-critical  component)  at  any  service  time 
(e.g.,  see  Fig.  1-3).  The  extent  of  damage  can  be  described 
quantitatively  by  the  number  of  structural  details  or  ligaments 
expected  to  exceed  specified  crack  size  limits  with  given  prob¬ 
ability  at  any  service  time  T  .  Hence,  the  durability  analysis 
provides  a  quantitative  description  of  structural  durability  in 
physically  meaningful  terms. 

Durability  analysis  procedures  are  summarized  as  fol¬ 
lows:  (1)  define  the  equivalent  initial  flaw  size  distribu¬ 

tion  (EIFSD)  using  fractographic  data  in  the  small  crack  size 
region  (e.g.,  0.01”  -  0.05"),  (2)  use  fractographic  data 

pooling  procedure  and  statistical  scaling  technique  to  esti¬ 
mate  the  EIFSD  parameters  in  a  "global  sense"  for  a  "single 
hole  population"  basis,  and  (3)  use  the  two-segment  determin¬ 
istic-stochastic  crack  growth  approach  (DCGA-SCGA)  [2,3,6]  to 
predict  the  extent  of  damage  in  the  entire  durability  criti¬ 
cal  component;  the  two-segment  deterministic  crack  growth 
approach  (DCGA-DCGA)  [2,3,  6,7]  is  also  reasonable  but  it  is 
slightly  less  conservative  than  the  DCGA-SCGA. 

Procedures  have  been  developed  for  defining  initial  fa- 
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J 

NOTES:  (l)  Largest  initial  flaw  size  in  each  structural  detail 
is  a  random  variable 

(2)  An  EIFSD  is  required  for  each  structural  detail  type 


Figure  1-2.  Initial  Fatigue  Quality  for  Each  Detail  Type 
is  Represented  by  an  Equivalent  Initial  Flaw 
Size  Distribution  (EIFSD). 
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AVE.  NO.  OF  DETAILS 
EXCEEDING  SPECIFIED 
CRACK  SIZE  LIMITS 


Figure  1- 


Functional  Impairment 


•  Extent  of  Damage  Concept  for  Quantitative 
Durability  Analysis. 


tigue  quality.  These  procedures  could  be  used  to  standardize 
the  way  initial  flaw  sizes  are  determined  from  fractographic 
data.  A  better  understanding  of  initial  flaw  sizes  (i.e., 
what  they  are  and  limitations)  has  been  developed  [e.g.,  2- 
4].  For  consistent  durability  analysis  predictions,  equiva¬ 
lent  initial  flaws  must  be  used  in  the  same  context  for  which 
they  were  defined.  This  means  that  equivalent  initial  flaws 
must  be  grown  forward  in  the  same  manner  the  EIFSs  were  es¬ 
tablished  by  back-extrapolating  fractographic  results. 
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SECTION  II 


DURABILITY  DESIGN  REQUIREMENTS 
AND  ANALYSIS  CRITERIA/GUIDELINES 

2 . 1  INTRODUCTION 

Tho  purpooe  <-ff  this  section  is  to:  (1)  briefly  review 
and  interpret  the  important  elements  of  the  Air  Force's  dur¬ 
ability  design  requirements  [8-11],  (2)  discuss  durability 

critical  parts  criteria,  (3)  provide  guidelines  and  recom¬ 
mended  formats  for  defining  quantitative  economic  life  cri¬ 
teria,  and  (4)  discuss  functional  impairment  due  to  fuel 
leaks  and  ligament  breakage. 

2.2  DURABILITY  DESIGN  REQUIREMENTS 
2.2.1  Objective  and  Scope 

The  objective  of  the  Air  Force  durability  design  re¬ 
quirements  [8-11]  is  to  minimize  in-service  maintenance  costs 
and  maximize  operational  readiness  through  proper  selection 
of  materials,  stress  levels,  design  details,  inspections,  and 
protection  systems.  These  design  requirements  include  both 
analyses  and  tests. 


2.2.2  General  Requirements 

Essential  durability  requirements,  conceptually  describ¬ 
ed  in  Fig.  2-1,  are  as  follows: 

o  The  economic  life  of  the  airframe  must  exceed  one 
design  service  life. 

o  No  functional  impairment  (e.g.,  loss  of  stiffness, 
loss  of  control  effectiveness,  loss  of  cabin  pres- 
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sura  or  fuel  leaks)  shall  occur  in  less  than  one 
design  service  life. 

o  The  economic  life  of  the  airframe  must  be  demon¬ 
strated  analytically  and  experimentally. 

2.2.3  Analytical  Requirements 

Analyses  are  required  to  demonstrate  that  the  economic 
life  of  the  airframe  is  greater  than  the  design  service  life 
when  subjected  to  the  design  service  loads  and  design  chemi¬ 
cal/thermal  environments.  The  economic  life  analysis  must 
account  for  initial  quality,  environment,  load  sequence,  ma¬ 
terial  property  variations,  etc.  The  analysis  must  be  veri¬ 
fied  by  tests. 


2.2.4  Experimental  Requirements 

Design  development  testa  are  required  to  provide  an 
early  evaluation  of  the  durability  of  critical  components  and 
assemblies  as  well  as  the  verification  of  the  durability  an¬ 
alysis. 

A  durability  test  of  a  full-scale  airframe  may  also  be 
required  by  the  Air  Force.  The  requirements  for  this  test 
are: 

1.  The  airframe  must  be  durability  tested  to  one  life¬ 
time.  Critical  structural  areas  must  be  inspected 
before  the  full  production  go-ahead  decision. 

2.  Two  lifetimes  of  durability  testing  plus  an  inspec¬ 
tion  of  critical  structural  areas  must  be  completed 
prior  to  delivery  of  the  first  production  aircraft. 

If  the  economic  life  of  the  airframe  is  not  reached  be- 
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fora  two  lifetimes  of  durability  testing,  the  following  op¬ 
tions  are  available: 

1.  Terminate  the  durability  testing  and  perform  a  non¬ 
destructive  inspection  followed  by  destructive 
teardown  inspection. 

2.  Terminate  the  durability  testing  and  perform  damage 
tolerance  testing  and  nondestructive  inspection 
followed  by  a  destructive  teardown  inspection. 

3.  Continue  the  durability  testing  for  an  approved 
period  of  time  followed  by  either  of  the  preceding 
options. 

2.3  DURABILITY  ANALYSIS  CRITERIA 
2.3.1  Durability  Damage  Modes 

There  are  several  modes  of  durability  damage,  including 
fatigue  cracking,  corrosion,  wear,  etc.  Due  to  its  impor¬ 
tance  and  prevalence,  fatigue  cracking  is  the  form  of  struc¬ 
tural  degradation  considered  in  this  handbook. 

2.3.2  Durability  Critical  Parts  Criteria 

Criteria  must  be  developed  for  determining  which  parts 
of  an  aircraft  are  durability  critical  (i.e.,  which  parts 
must  be  designed  to  meet  the  durability  design  requirements) . 
The  durability  critical  parts  criteria  vary  from  aircraft  to 
aircraft.  They  are  especially  dependent  on  the  definition  of 
economic  life  for  the  particular  aircraft  involved.  A  typi¬ 
cal  flow  diagram  for  selecting  which  parts  are  durability 
critical  is  presented  in  Fig.  2-2.  In  Fig.  2-2,  durability 
refers  to  the  ability  of  an  airframe  to  resist  cracking 
whereas  damage  tolerance  refers  to  the  ability  of  an  airframe 


2-4 


loo  STATIC  STR 

s 

i 

1 

O ART  ALWAYS  LOAOf 

SIM  CQMPRCSIION?> 

VI* 


|00  PURAHLITY  AMALVtttl  ICO  0AMA01 TOLIHAMCIA  OURAOIUTV  ANAIySs] 


1 


JL 


Cpumbuty  COWTROI  meg  livut) 

1v5  TSo 


< 


6aMA8K  TQIIRANCS  MOM  CRITICAL  V 

THAW  QUMKlim  / 


r 

dwa  cowTwoLiio  MATiRUift 
^JMO  ' 

db 


<; 


E. 


0AMA8C  TQLIRAMCI  COWTMlN^ 
STRItt  LKVIL? 


I 


MO 


dm  controllip  matiwai^ 
1*q  — 


|UM  COMTBOILI 0  MATIRIAtl 


nnnrw 


MO 


c 


US<  CONTROLLED  MATIMAL 


yilNOLl  LOAO  PATH  OR 
"\JIOW|W*eTMUT 


ado  to 


1 


tit 


> 


ADO  TO  FRACTURS  CRITICAL 
PARTI  LIST 


fllih 


Figure  2-2.  Flow  Diagram  for  Selecting  Durability  Critical 
*  Parts. 
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to  resist  failure  due  to  the  presence  of  such  cracks. 


2.3.3  Economic  Life  Criteria/ Guidelines 

Criteria  must  be  developed  for  determining  the  economic 
life  of  the  particular  aircraft  of  interest.  Similar  to  the 
durability  critical  parts  criteria,  economic  life  criteria 
vary  from  aircraft  to  aircraft.  They  may  be  based  on  fasten¬ 
er  hole  repair  (e.g.,  reaming  the  damaged  fastener  hole  to 
the  next  nominal  hole  size),  functional  impairment  (e.g., 
fuel  leakage),  residual  strength,  etc.  Two  promising  analy¬ 
tical  formats  for  quantifying  the  economical  life  of  an  air¬ 
frame  are  (1)  the  probability  of  crack  exceedance,  and  (2) 
cost  ratio:  repair  cost/replacement  cost.  Both  formats  re¬ 
quire  a  durability  analysis  methodology  capable  of  quantify¬ 
ing  the  extent  of  aircraft  structural  damage  as  a  function  of 
service  time.  For  example,  assume  the  economic  life  criteria 
are  based  on  the  number  of  fastener  holes  which  cannot  be 
economically  repaired  (i.e. ,  number  of  fastener  holes  with 
crack  sizes  equal  to  or  greater  than  specified  size  xx) . 
Then  an  analytical  format  for  quantifying  economic  life  is 
presented  in  Fig.  2-3.  In  Fig.  2-3,  P  is  the  exceedance 
probability.  Various  aspects  of  economic  life  are  discussed 
further  in  the  following  subsections  and  elsewhere  [12-23], 

2.3.3. 1  Economic  Life  Definition 

The  economic  life  of  an  aircraft  structure  is  currently 
defined  in  qualitative  terms:  "...the  occurrence  of  wide¬ 
spread  damage  which  is  uneconomical  to  repair  and,  if  not  re¬ 
paired,  could  cause  functional  problems  affecting  operational 
readiness"  [8,10].  Acceptable  limits  for  "widespread  damage" 
and  "uneconomical  repairs"  must  be  defined  for  each  aircraft 
design  and  such  limits  must  be  approved  by  the  Air  Force. 

A  quantitative  definition  of  economic  life  is  not  given 
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Figure  2-3.  Analytical  Format  for  Economic  Life. 
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in  this  handbook.  However,  guidelines  are  presented  for 
specifying  economic  life  criterion  (Ref.  Section  2. 3. 3. 4). 
In  any  case,  quantitative  criteria  for  the  economic  life  of 
aircraft  structures  should  be  based  on  specific  aircraft  re¬ 
quirements  and  the  user's  acceptable  limits  for  aircraft  per¬ 
formance  and  maintenance  costs. 


2. 3. 3. 2  Economic  Repair  Limit 

The  "economic  repair  limit"  is  the  maximum  crack  size  in 
a  structural  detail  that  can  be  economically  repaired.  Such 
limits  can  easily  be  defined  from  geometric  considerations 
for  fastener  holes  but  such  limits  are  more  difficult  to  de¬ 
fine  for  structural  details  such  as  cutouts,  fillets,  etc. 
For  example,  the  economic  repair  limit  for  a  fastener  hole 
may  be  governed  by  the  largest  radial  crack  that  can  be 
cleaned  up  by  reaming  the  hole  to  the  next  fastener  size 
(e.g.,  0.03"  to  0.05"  radial  crack). 

The  objective  of  the  durability  analysis  method  present¬ 
ed  in  this  handbook  is  to  analytically  predict  the  number  of 
structural  details  with  a  crack  size  which  would  cause  an  un¬ 
economical  repair  or  functional  impairment.  The  user  must 
define  the  uneconomical  repair  or  functional  impairment  crack 
size  for  the  details  to  be  included  in  the  extent-of -damage 
assessment.  Such  crack  sizes  depend  on  considerations  such 
as  structural  detail  type,  location,  accessability,  inspect- 
ability,  repairability,  repair  costs,  etc. 

Structural  details  may  contain  one  or  more  cracks.  How¬ 
ever,  structural  durability  is  concerned  with  the  largest 
crack  in  each  detail  which  may  require  repair  or  part  re¬ 
placement  . 


2 . 3 . 3 . 3  Extent  of  Damage 


The  extent  of  damage  is  a  quantitative  measure  of  the 
number  of  structural  details  containing  cracks  that  exceed 
specified  crack  size  limits  as  a  function  of  service  time. 
Structural  maintenance  requirements  and  costs  depend  on  the 
number  of  structural  details  requiring  repair.  The  '•durab¬ 
ility"  of  the  structure  depends  on  the  extent  of  damage  for 
the  population  of  structural  details  in  a  part,  a  component, 
or  airframe. 

The  statistics  of  the  extent  of  damage,  such  as  mean  and 
extreme  values  for  selected  probabilities,  can  be  predicted 
using  the  analytical  tools  provided  in  this  handbook.  Extent 
of  damage  predictions  provide  the  basis  for  analytically  en¬ 
suring  that  a  durability-critical  part  or  component  will  not 
crack  excessively  in  less  than  one  service  life. 

2. 3. 3. 4  Formats  for  Economic  Life  Criteria 

Two  analytical  formats  for  defining  quantitative  econom¬ 
ic  life  criteria  are  recommended:  (1)  probability  of  crack 
exceedance^  anc  (2)  cost  ratio:  repair  cost/replacement  cost 
[15-18].  The  analytical  tools  described  in  this  handbook  can 
be  used  to  predict  results  in  these  formats.  Various  aspects 
of  each  format  for  a  quantitative  economic  life  criterion  are 
discussed  below,  including  examples  and  guidelines  (Ref.  Fig. 
2-3)  . 

2. 3. 3. 4.1  Probability  of  Crack  Exceedance.  The  probability 
that  a  crack  will  be  larger  than  a  specified  crack  size  at  a 
particular  service  time  is  referred  to  as  the  ’’probability  of 
crack  exceedance."  This  quantity  is  a  fundamental  output  of 
the  durability  analysis  methodology  described  in  this  handbook. 
For  example,  in  Fig.  2-4  the  probability  of  exceeding  crack 
size  X]_  at  t  =  T*  is  represented  by  the  cross-hatched  area  un- 
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der  the  crack  size  density  function  at  t  -  T  .  When  the  de¬ 
terministic  crack  growth  approach  is  used,  crack  size  rank¬ 
ings  in  the  respective  distributions  for  two  different  times 
are  preserved;  namely,  the  crack  size  x^  at  t  *  T  has  the  same 
rank  (or  percentile)  as  the  initial  crack  size  at  y^(T*)  at  t 
*  0.  The  probability  of  crack  exceedance  can  be  used  to  pre¬ 
dict  the  expected  number  of  repairs  in  a  given  service  interval 
[16,18].  It  also  provides  a  basis  for  judging  airframe  dur¬ 
ability  and  for  analytically  demonstrating  design  compliance 
with  the  governing  criterion  for  economic  life. 


Another  explanation  of  the  probability  of  crack  exceed¬ 
ance  concept  will  now  be  given.  Each  common  structural  de¬ 
tail,  in  a  group  of  details  having  a  common  stress  history, 
has  a  single  dominant  crack.  Such  cracks  form  a  crack  popu¬ 
lation  and  their  "initial"  size  depends  on  the  manufacturing 
quality  for  each  structural  detail.  The  probability  of  exc¬ 
eeding  crack  size  x^  at  time  T  is  represented  by  the  cross- 
hatched  area  under  the  probability  density  of  crack  sizes 
shown  in  Fig.  2-4.  Suppose  the  probability  of  crack  exceed¬ 
ance  is  p(i,‘f')  -  0.05.  This  means  that  on  the  average  5%  of 
the  details  (e.g.,  5%  of  the  fastener  holes)  in  a  part  or 
component  would  be  expected  to  have  a  crack  size  >  x^  at 
time  T  .  P(i,7~)  is  a  fundamental  measure  of  the  extent  of 
damage.  Using  the  binomial  distribution,  the  extent  of  dam¬ 
age  for  different  groups  of  details  can  be  combined  to  quan¬ 
tify  the  overall  damage  for  a  part,  a  component  or  airframe. 


The  allowable  crack  exceedance  is  one  criterion  recom¬ 
mended  for  quantifying  economic  life.  Although  this  handbook 
provides  guidelines  for  quantifying  the  allowable  crack  ex¬ 
ceedance,  specific  values  are  not  presented  for  demonstrating 
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design  compliance  with  the  Air  Force's  durability  design  re¬ 
quirements.  Such  values  must  be  tailored  for  specific  air¬ 
craft  structure  and  the  user's  acceptable  limit  for  structur¬ 
al  maintenance  requirements/costs,  functional  impairment,  op¬ 
erational  readiness,  etc.  The  allowable  crack  exceedance 
criterion  for  economic  life  design  compliance  shall  be  ap¬ 
proved  by  the  Air  Force. 

The  allowable  crack  exceedance  for  a  part  or  component 
depends  on  several  factors,  including:  criticality,  access- 
ability,  inspectability,  repairability,  cost,  operational 
readiness,  acceptable  risk  limits,  etc.  For  example,  an  ex¬ 
pensive  fracture  critical  part  may  be  embedded  into  the  wing 
under-structure.  The  part  is  not  readily  accessible  and  it 
is  difficult  to  inspect  and  repair.  Suppose  the  bolt  hole 
for  this  part  governs  its  economic  life.  Then  a  lower  allow¬ 
able  crack  exceedance  may  be  desired  than  for  an  equally  cri¬ 
tical  part  that  is  more  accessible  and  inspectable.  For  ex¬ 
ample,  an  average  of  2t  crack  exceedance  at  1.2  service  lives 
might  be  suitable  in  the  first  case  and  an  average  of  5% 
might  be  appropriate  for  different  circumstances. 

An  example  for  the  probability  of  crack  exceedance  cri¬ 
terion  is  as  follows.  The  economic  life  of  a  part  or  compo¬ 
nent  is  reached  when  5  percent  of  the  structural  details 
(e.g.,  fastener  holes,  cutouts,  fillets,  etc.)  have  reached  a 
crack  size  >  a  specified  limiting  crack  size  at  1.2  service 
lives.  The  limiting  crack  size  depends  on  the  type  of  struc¬ 
tural  detail,  the  economic  repair  limit,  and  the  crack  size 
which  would  cause  functional  impairment  (limiting  case) . 
Structural  safety  or  damage  tolerance  must  not  be  compromis¬ 
ed.  Also,  the  specified  limiting  crack  size  for  each  detail 
type  should  account  for  inspection  capabilities  and  require¬ 
ments,  and  operational  readiness . 

The  economic  life  criterion  described  (i.e.,  5%  crack 
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exceedance)  can  ba  used  to  demonstrate  economic  life  design 
compliance  analytically  and  experimentally.  The  analytical 
tools  presented  in  this  handbook  can  be  used  to  quantify  the 
extent  of  damage  in  terms  of  crack  exceedance.  Therefore, 
given  the  criterion  for  economic  life,  design  compliance  can 
be  analytically  assured.  Experimental  compliance  can  be  de¬ 
termined  based  on  the  results  of  the  durability  demonstration 
test  results. 


2 . 3 . 3 . 4 . 2  Repair  Cost/Replacement  Cost  Ratio.  The  ra¬ 
tio  of  repair  cost/replacement  coat  is  another  recommended 
criterion  for  quantitative  economic  life.  For  example,  when 
the  cost  to  repair  a  part  or  component  exceeds  the  cost  to 
replace  it,  the  economic  life  is  reached.  In  other  words, 
the  economic  life  is  reached  when  the  cost  ratio  -  1  at  a 
specified  service  life  (e.g.,  1.2  service  lives). 

Input  from  the  aircraft  user  is  needed  to  define  accept¬ 
able  allowable  cost  ratios  for  different  parts  or  components. 
Allowable  cost  ratios  could  be  specified  for  particular  de¬ 
sign  situations  and  user's  goals. 

Repair  costs  are  proportional  to  the  number  of  structur¬ 
al  details  (e.g.,  fastener  holes)  requiring  repair  after  a 
specified  service  time.  The  analytical  tools  described  in 
this  handbook  can  be  used  to  quantify  the  number  of  details 
requiring  repair  as  a  function  of  service  time.  Although 
specific  repair  cost  data  may  be  difficult  to  obtain  for  dif¬ 
ferent  circumstances  and  replacement  costs  may  vary,  the  cost 
ratio  can  be  estimated  using  assumed  repair  and  replacement 
costs. 

The  cost  ratio  criterion  for  economic  life  is  not  recom¬ 
mended  for  demonstrating  design  compliance  unless  acceptable 
cost  data  are  available.  However,  this  criterion  is  recom¬ 
mended  for  evaluating  user  design  tradeoff  options  affecting 


2-13 


the  life-cycle-cost  of  the  airframe.  The  analytical  tools 
dssoribad  in  this  handbook  can  bs  ussd  to  svaluats  ths  life- 
cycle-cost  design  tradeoffs. 

2.4  FUNCTIONAL  IMPAIRMENT  DUE  TO  FUEL  LEAKS/ LIGAMENT 

BREAKAGE 

Fuel  leaks  and  ligament  breakage  are  other  forms  of 
functional  impairment  which  must  be  accounted  for  in  the  de¬ 
sign  of  metallic  airframes.  Large  through-the-thickness 
cracks  may  cause  fuel  leaks  with  a  progressive  increase  in 
the  state  of  damage  as  a  function  of  service  time.  Such 
cracks  may  not  pose  an  immediate  structural  safety  problem. 
However,  they  may  affect  the  operational  readiness  of  the 
aircraft  and  increase  the  structural  maintenance  requirements 
and  repair  costs.  Without  inspection,  repair  and  maintenance 
the  damaged  areas  may  eventually  lead  to  a  safety  of  flight 
problem  or  necessitate  expensive  repairs, 

Fuel  leaks  are  a  fire  hazard.  They  increase  fuel  con¬ 
sumption  and  may  affect  the  operational  readiness  of  the  air¬ 
craft.  Ligament  breakage  is  like  a  cancer  in  that  damage  may 
continue  to  spread  with  increased  service  time.  For  example, 
a  crack  in  one  hole  may  grow  to  an  adjacent  hole  and  may  con¬ 
tinue  to  grow  in  service  to  other  holes.  An  adjacent  hole 
may  act  as  a  crack  stopper  but  the  redistribution  or  shedding 
of  load  around  the  damaged  area  under  service  conditions  may 
eventually  spread  the  damage  to  adjacent  areas. 

The  durability  analysis  tools  in  this  handbook  can  be 
used  to  assess  the  probability  of  functional  impairment  due 
to  excessive  cracking,  fuel  leaks  or  ligament  breakage.  Such 
tools  can  be  used  during  the  design  stage  to  select  mate¬ 
rials,  design  concepts  and  allowable  stress  levels  to  satisfy 
the  Air  Force's  durability  design  requirements. 
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SECTION  III 

SUMMARY  OF  THE  DURABILITY  ANALYSIS  METHOD 

Essential  elements  and  aquations  of  tha  durability  an¬ 
alysis  msthod  ara  summarized  in  this  section.  The  technical 
approach,  based  on  the  two-segment  DCGA-SCGA,  is  presented, 
including  stap-by-step  procedures  for  implementing  the  me¬ 
thod.  Durability  analysis  guidelines  are  given  in  Section  IV 
and  the  methods  are  demonstrated  in  Section  V.  Further  de¬ 
tails  are  given  elsewhere  [2-7,24-29]. 

3.1  GENERAL  DESCRIPTION  OF  THE  TECHNICAL  APPROACH 

The  initial  fatigue  quality  (IFQ)  of  structural  details 
reflected  in  tha  durability  analysis  is  represented  by  an 
equivalent  initial  flaw  size  distribution  (EIFSD) .  An  EIFSD 
is  defined  by  structural  detail  types  (e.g.,  see  Fig.  1-2). 
An  equivalent  initial  flaw  size  (E1FS)  is  an  artificial  crack 
size  at  time  zero  which  results  in  an  actual  crack  size  when 
the  EIFS  is  grown  forward.  EIFSs  are  determined  by  back-ex¬ 
trapolating  fractographic  results  to  time  zero.  An  EIFS  has 
no  direct  relationship  to  actual  initial  flaws  (e.g., 
scratches,  burrs,  microdefects,  etc.)  in  a  structural  detail 
and  such  flaws  cannot  be  verified  by  NDI.  An  EIFS  is  assumed 
to  be  a  random  variable  which  is  statistically  described  by 
an  EIFSD. 

Once  a  suitable  EIFSD  has  been  determined,  the  EIFSD  is 
grown  forward  using  a  two-segment  deterministic-stochastic 
crack  growth  rate  model.  The  probability  of  crack  exceed¬ 
ance,  p(i,0,  at  any  service  time  T  and/or  the  cumulative 
distribution  of  service  time  to  reach  any  specified  crack 
size  x1#  Ft^x  j  (T) ,  can  be  predicted  for  a  durability-criti¬ 
cal  component.1  Using  the  probability  of  crack  exceedance 
predictions,  the  assumption  of  statistically  independent 
cracking  and  the  binomial  distribution,  the  extent  of  damage 
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mean  (P  ■  0.5)  and  uppar  bound  limit  (e.g.,  P  -  0.05)  can  ba 
aatiaatad  for  salactad  exceedance  probabilitias  F.  The  ex¬ 
tant  of  damage  dafinaa  statistically  the  number  of  structural 
details  or  ligaments  in  a  durability-critical  component  ex¬ 
pected  to  exceed  crack  size  limits  for  functional  impairment 
at  a  given  service  time.  Hence,  the  extent  of  damage  pro¬ 
vides  a  basis  for  incorporating  durability  requirements  into 
the  design  process  and  for  evaluating  durability  design 
tradeoffs  (e.g.,  material,  design  concept,  stress  level,  load 
spectra,  %  bolt  load  transfer,  etc.). 

The  technical  approach  for  the  durability  analysis  in¬ 
cludes  four  essential  steps:  (1)  determine  the  initial  fa¬ 
tigue  quality  or  EIFSD  suitable  for  the  structural  details  to 
be  reflected  in  the  durability  analysis,  (2)  determine  a 
suitable  service  crack  growth  master  curve  in  two  crack 
growth  segments,  (3)  predict  the  probability  of  crack  exceed¬ 
ance  p(i, T)r  et  a  given  service  time,  and/or  the  cumulative 
distribution  of  service  time  to  reach  a  given  crack  size  x^ 
and  (4)  estimate  the  extent  of  damage  mean  (P  *  0.5)  and  up¬ 
per  bound  limit  for  selected  exceedance  probability  (e.g.,  P 
■  0.05)  at  any  given  service  time. 

Essential  elements  and  features  of  the  durability  analy¬ 
sis  method  are  conceptually  described  in  Figs.  3-1  through  3- 
5.  Details  are  given  in  the  following  and  elsewhere  [2-7, 
24-29] . 

The  durability  analysis  method,  based  on  the  two-segment 
DC6A-5CGA,  has  been  demonstrated  for  clearance-fit  fastener 
holes  in  7475-T7351  aluminum  for  both  protruding  head  and 
countersunk  fasteners.  The  method  has  been  demonstrated  for 
both  small  (e.g.,  <  0.05")  and  large  through-the-thickness 
(e.g.,  0.5"  -  0.75")  fatigue  cracks.  The  durability  analysis 
method  presented  needs  to  be  further  investigated  for  other 
structural  details  (e.g.,  lugs,  cutouts,  fillets,  etc.)  and 
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Element 3  of  the  Durability  Analysis  Mathodc-logy. 


(a)  Probability  of  Crack  Exceedance  at  Service  Tima  T 


(b)  Cumulative  Distribution  of  Sarvica  Tima  at 
Crack  Size  x1 . 


Figure  3-3*  Qoodnass-Of-Flt  Plots  for  Justifying  an 
EIFSD  for  Durability  Analysis. 


3-5 


fo(t)/dt  : 

QlU(t)Jbl 

SEMIICE  HIE 


Figure  3-4.  Two-Segment  DCGA-SCGA  for  Durability 
Analysis. 
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P  :  Exceedance  Probability 


Figure  3-5.  Extent  of  Damage  Mean  and  Upper  Bound 
Limit  Concept. 
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for  fastener  holes  with  fatigue  life  enhancement  such  as  in¬ 
terference  fit  fasteners,  cold-working,  etc.  In  any  case,  a 
general  method  and  basic  framework  for  performing  ''quantita¬ 
tive"  durability  analysis  has  been  established. 

3.2  INITIAL  FATIGUE  QUALITY 

The  initial  fatigue  quality  <IFQ)  defines  the  initially 
manufactured  state  of  a  structural  detail  or  details  with 
respect  to  initial  flaws  in  a  part,  component,  or  airframe 
prior  to  service.  The  IFQ  for  a  group  of  replicate  details 
(e.g.,  fastener  holes)  is  represented  by  an  equivalent  in¬ 
itial  flaw  size  (EIFS)  distribution. 

The  Weibull  compatible  distribution  function  proposed  by 
Yang  and  Manning  [16,30]  has  been  found  to  be  reasonable  for 
representing  the  EIFS  cumulative  distribution  [1,15-22,25-30] 

(3-1) 

«  /.a  j  X  > 

in  which  xu  -  EIFS  upper  bound  limit;  cC  and  f  are  empiri¬ 
cal  parameters. 

An  EIFS  value  for  a  fastener  hole  is  determined  by  back- 
extrapolating  fractographic  data  in  a  selected  crack  size 
range  (AL-AU)  using  a  simple  but  versatile  deterministic 
crack  growth  rate  model  recommended  by  Yang  and  Manning  [16, 
30], 

y 

<?[«.<*>]  (J_2) 

where  da(t)/dt  -  crack  growth  rate,  a(t)  -  crack  size  at  any 
time  t  in  flight  hours,  and  Q  and  b  are  empirical  crack 
growth  rate  parameters.  The  special  case  b  ■  1  is  used  here¬ 
in. 
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After  EIFS  values,  a(0) ,  are  obtained  from  all  available 
fractographic  data,  they  are  fitted  by  Eq.  (3-D  to  determine 
the  EIFS  distribution  (EIFSD)  parameters  xu,  Pi,  and  p  .  To 
predict  the  extent  of  cracking  in  service,  the  equivalent  in¬ 
itial  flaw  size  distribution  is  grown  forward,  and  the  dis¬ 
tribution  of  the  crack  size  a(t)  at  any  service  time  t  can  be 
derived  from  that  of  a(0)  given  by  Eq.  C3-1).  The  EIFSD  is 
grown  forward  to  predict:  (1)  the  probability  that  a  crack 
in  the  ith  stress  region  at  any  service  time,  7* ,  will  exceed 
any  given  crack  size,  x1#  denoted  by  p(i,T),  and  (2)  the  cum¬ 
ulative  distribution  of  service  time,  j  (T)  #  for  a  crack 

in  the  ith  stress  region  to  reach  any  given1  crack  size 
p(i,  T)  is  referred  to  as  the  crack  exceedance  probability. 
The  two-segment  deterministic-stochastic  crack  growth 
approach  (DCGA-SCGA)  is  described  in  Section  3.4  and  in  Fig. 
3-4  for  growing  the  EIFSD  forward  to  predict  p(i,T)  and/or 


3.3  ESTIMATION  AND  OPTIMIZATION  OF  EIFSD  PARAMETERS 

The  essential  procedures,  concepts,  and  equations  for 
estimating  and  optimizing  EIFSD  parameters  xu, oCand  p  in  Eq. 
3-1  are  given  in  the  following.  Details  are  given  elsewhere 
[2,3].  Six  major  topics  are  covered:  (1)  general  procedure, 
(2)  determination  of  EIFSs,  (3)  data  pooling,  (4)  statistical 
scaling  technique,  (5)  combined  least  square  sums  approach 
(CLSSA)  and  (6)  optimization  of  parameters  and  goodness-of- 
-fit. 

3.3.1  General  Procedure 

1.  Select  fractographic  data  set(s)  to  be  used  to  de¬ 
termine  the  EIFSD.  The  data  sets  should  be  for  the  same 
material,  same  type  load  spectrum  (e.g.,  fighter,  bomber  or 
transport)  and  type  fastener/hole/fit  (i.e.,  straight  bore  or 
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countersunk) . 


2.  Screen  each  fractographic  data  set  for  "anomalies" 
using  software  filename  =  "SCREEN"  from  Volume  V  [24],  Cen¬ 
sor  out  fatigue  cracks  with  anomalous  behavior  (e.g.,  surface 
crack  rather  than  crack  in  bore  of  hole,  cracks  with  crack 
growth  rate  extremes  (i.e.,  very  fast  or  very  slow  compared 
to  rest  of  data,  etc.)* 

3.  Select  a  suitable  fractographic  crack  size  range, 
AL-AU  (e.g.,  .01"-. 05"),  and  a  reference  crack  size,  x^,  for 
defining  service  times  or  time-to-crack  initiation  (TTCI) . 

4.  Estimate  the  EIFSs  for  each  screened  data  set  as 
described  in  Section  3.3.2.  Use  the  largest  fatigue  crack  in 
each  specimen  to  determine  the  EIFSs. 

l/  5.  Assume  the  EIFSD  is  represented  by  the  Weibull  com¬ 
patible  distribution  function  given  in  Eq.  (3-1).  Other  EIFSD 
functions  could  also  be  used  if  appropriate  (e.g.,  lognormal 
compatible  [2],  lognormal,  two  parameter  Weibull,  etc.).  For 
a  given  EIFS  upper  bound  limit,  xu  (e.g.,  largest  EIFS  in 
data  set(s)  <  xu  <  0.05")  estimate  the  EIFSD  parameters  od  and 
in  Eq.  0-1)  using:  (1)  the  EIFSs  from  Step  4,  (2)  data  pool¬ 
ing  (Section  3.3.3),  (3)  statistical  scaling  (Section  3.3.4), 
and  (4)  combined  least  square  sums  approach  (CLSSA)  and  an 
"EIFS  fit"  (Section  3.3.5). 

6.  Optimize  the  EIFSD  parameters  using  the  iterative 
procedure  described  in  Section  3.3.6.  For  clearance-fit 
fastener  holes  and  the  Weibull  compatible  distribution  func¬ 
tion,  an  EIFS  upper  bound  limit  range  of  xu  =  0.02"-0.05"  is 
reasonable. 

7.  Justify  the  candidate  EIFSD  for  durability  analysis 
by  checking  goodness-of-fit  (see  Fig.  3-3) .  Correlate  theo- 
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retical  predictions  for  the  cumulative  distribution  of  ser¬ 
vice  time  to  reach  the  crack  size  ^  (T) ,  and/or  the 
cumulative  distribution  of  crack  size  at  service  time  T, 

Y m  (x) ,  with  fractographic  results  for  selected  data  sets. 
Check  goodness-of-fit  using  those  fractographic  data  sets 
that  were  used  to  estimate  the  EIFSD.  Where  possible/  other 
fractographic  data  sets  should  also  be  used  to  check  the 
"goodness-of-fit . " 

Computer  software,  briefly  described  in  Section  VI,  is 
available  in  Volume  V  [24]  to  estimate,  optimize  and  justify 
the  EIFSD  parameters  for  durability  analysis,  including  a 
goodness-of-fit  plotting  capability.  This  software  can  be 
implemented  on  an  IBM  or  IBM-compatible  personal  computer. 

3.3.2  Determination  of  EIFSs 

Equivalent  initial  flaw  sizes  (EIFSs)  for  a  given  frac¬ 
tographic  data  set  are  determined  by  back-extrapolating  frac¬ 
tographic  results.  The  EIFS  master  curve  for  each  fracto¬ 
graphic  data  set  is  defined  by  integrating  Eq. (3-2) (with  b«l) 
from  a(0)  to  a(t)  to  obtain 

a(t)  «  a(0)  exp(Qt)  (3-3) 


or 


EIFS  -  a(0)  -  a(t)  exp(-Qt) 


(3-4) 


in  which  Q  -  empirical  crack  growth  rate  parameter  (referred 
to  as  "pooled  Q"  for  a  data  set),  a(t)  -  crack  size  at  any 
time  t,  a(Q)  -  EIFS  »  crack  size  at  t-0. 


3-11 


The  "pooled  "Q"  value  in  Eq.  (3-3)  or  (3-4)  for  a  data  set 
can  be  determined  as  follows.  Suppose  the  ith  fractographic 
data  set  contains  a  total  of  m  fatigue  cracks,  where  each  fa¬ 
tigue  crack  is  denoted  by  j  -  1,  2,  . m.  The  jth  fatigue 
crack  has  a  total  of  N  pairs  of  fractographic  data  in  the  AL- 
AU  range,  denoted  by  [^(t^),*  t^] ,  i.e.,  a^t^)  -  kth  crack 
size  for  the  jth  fatigue  crack  at  service  time  t^  in  the  AL- 
AU  range,  where  k  «  l,  2,  . ..,  N. 


The  crack  growth  rate  parameter  for  a  single  fatigue 
crack,  say  the  jth  fatigue  crack,  denoted  by  Qj ,  is  estimated 
from  fractographic  data  of  the  jth  fatigue  crack  in  the  AL-AU 
range  using  Eq.  ,(3-3)and  the  least  squares  fit  procedure  as 
follows 


(3-5) 


in  which  Xk  -  tk  and  Y^  -  In  aj (t^) . 

Qj  given  in  Eq.  (3-5)  is  the  crack  growth  rate  parameter 
for  the  jth  crack  and  it  is  obtained  using  the  fractographic 
data  of  the  jth  crack.  Let  be  the  crack  growth  rate  para¬ 
meter  for  the  ith  data  set  consisting  of  *n  cracks.  Then, 
is  referred  to  as  the  "pooled  Q"  value  for  the  ith  data  set. 
It  is  obtained  using  all  the  fractographic  data  in  the  ith 
data  set,  i.e.,  all  fractographic  data  for  w  cracks  in  the 
AL-AU  range,  and  the  least  squares  fit  procedure, 


(3-6) 
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One*  C  has  bssn  determined,  ths  EIFSs  for  ths  ith  data 
sst  can  be  determined  as  follows.  Select  a  reference  crack 
size  x1  within  the  AL-AU  range  (i.e.,  AL  <  aQ  <  AU)  for  the 
time-to-crack -  initiation.  Then,  determine  for  each  fatigue 
crack  in  the  data  sat,  by  interpolation  or  extrapolation,  the 
time-to-crack-initiation  (TTCI)  and  the  results  are  denoted 
by  (Tr  T2„  ...,  Ta). 

The  EIFS  sample  value  for  the  jth  crack,  denoted  by 
aj(0),  is  obtained  from  Eq.  (.3-4) by  setting  t  s  Tj  and  a(t)  « 
a(Tj )  »  aQ,  i.e., 

aj(0)  -  jth  EIFS  -  aQ  exp(->  Q^Tj )  ;  j  -  1,  2,  ...,  m  (3-7) 

In  other  words,  the  jth  crack  with  an  EIFS,  aj(0),  will  grow 
to  the  reference  crack  size  aQ  at  Tj .  Equation  (3-7)  can  also 
be  interpreted  in  another  form;  i.e.,  the  reference  crack  size 
aQ  at  Tj  is  back-extrapolated  to  time  zero  to  determine  the 
corresponding  EIFS  value,  aj(0). 

Note  that  all  EIFS  values  for  the  ith  data  set,  i.e., 
[a1(0) ,a2(0) ,  ...,  aa(0)],  are  computed  from  Eq. (3-7)  using 
the  same  value;  i.e.,  the  "pooled  Q"  value  for  the  ith 
data  set.  The  EIFS  values  thus  obtained  are  referred  to  as 
"deterministic-based  EIFSs"  [2,3,6,7,27,29].  When  the  jth 
EIFS  value-  for  the  jth  crack  is  computed  using  its  own  crack 
growth  rate  parameter  Qj ,  Eq.  (3-5,)  the  EIFS  data  set  thus  es¬ 
tablished  is  referred  to  as  the  "stochastic-based  EIFS"  [2,3, 
5,25,27].  As  a  result  of  extensive  investigations  conducted 
[2,3,5-7,25,27,29,59],  the  deterministic-based  EIFS  has  been 
recommended. 
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Finally#  given  the  crack  growth  rate  parameter  or 
"pooled  QM  value  for  the  ith  data  set,  the  crack  size-time 
relationship  given  by  Eq.  (3-4), i.e.,  a(0)  -  aCtJexpC-Q^) ,  is 
referred  to  as  the  "EIFS  master  curve"  for  the  ith  data  set. 
In  fact,  the  back-extrapolation  for  computing  all  EIFS  sample 
values  for  the  ith  data  set  uses  the  same  EIFS  master  curve. 

3.3.3  Data  Pooling 

Suppose  we  have  M  different  fractographic  data  sets  gen¬ 
erated  under  different  test  conditions  (e.g.,  same  type  of 
load  spectrum  but  different  stress  levels  and  t  bolt  load 
transfer).  Then,  we  have  M  different  "pooled  QM  values; 
i.e.,  Qi  (i  -  1,  2,  ...,  M),  and  M  different  EIFS  master 

curves.  Consequently,  M  EIFS  data  sets  can  be  computed  using 
the  corresponding  "pooled  Q"  value  for  each  data  set  as  des¬ 
cribed  previously.  These  M  data  sets  of  EIFS  values  can  be 
pooled  together,  referred  to  as  "data  pooling"  (see  Fig.  3- 
2),  to  define  the  EIFSD  parameters.  Pooling,  effectively  in¬ 
creases  the  sample  size  and  confidence  in  the  EIFSD  paramet¬ 
ers.  Likewise,  this  is  a  'reasonable  approach  for  justifying 
an  EIFSD  for  mors  general  applications. 

The  fractographic  data  pooling  concept  is  conceptually 
illustrated  in  Fig.  3-6  using  two  fractographic  data  sets. 
The  premise  of  data  pooling  is  that  each  data  set  has  a  com¬ 
mon  EIFSD.  For  example,  if  the  TTCIs  for  each  data  set  are 
regressed  backwards  to  time  zero  using  the  applicable  EIFS 
master  curve  for  each  data  set,  the  resulting  EIFSs  have  the 
same  EIFSD.  In  other  words,  the  TTCI  distribution  for  dif¬ 
ferent  fractographic  data  sets  can  be  determined  using  the 
same  EIFSD. 

Thus,  EIFS  values  for  M  data  sets  obtained  previously 
can  be  pooled  together  and  used  to  determine  the  EIFSD  param¬ 
eters,  xu,  o i  and  ^  .  The  combined  least  square  sums  proce- 
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dure  to  be  described  later  nay  be  used  to  optimize  the  EXFSD 
parameters.  The  approach  described  above  is  referred  to  as 
the  "EIFS  fit"  [2,3],  Another  approach  to  determine  the 
EXFSD  parameters  is  a  direct  application  of  TTCI  data  sets  in 
conjunction  with  applicable  EIFS  master  curves  (or  back-ex¬ 
trapolation  transformation)  as  described  in  Refs.  2  and  3. 
Such  an  approach  is  referred  to  as  the  "TTCI  fit"  [2,3].  Al¬ 
though  there  are  subtle  differences  between  the  two  approach¬ 
es,  either  approach  gives  the  same  EXFSD  parameters  [3].  The 
two  approaches  are  conceptually  described  in  Fig.  3-7.  An 
"EIFS  fit"  is  recommended  and  is  emphasized  in  the  following 
because  EIFS  statistics  (i.e.,  mean  and  standard  deviation) 
provide  a  common  baseline  for  comparing  and  cataloging  "init¬ 
ial  fatigue  quality"  data  from  various  sources;  whereas,  TTCI 
statistics  (i.e.,  mean  and  standard  deviation)  do  not. 

Software  is  available  in  Volume  V  [24]  for  implementing 
the  data  pooling  procedure  described  in  this  section  on  an 
IBM  or  IBM-compatible  PC. 

3.3.4  Statistical  Scaling  Technique 

The  IFQ  or  EIFSD  for  fastener  holes  is  defined  for  a 
"single  hole  population."  Therefore,  the  fatigue  cracking 
resistance  of  each  fastener  hole  in  each  test  specimen  is  ac¬ 
counted  for  in  the  definition  of  the  EXFSD.  Test  specimens 
for  acquiring  fatigue  crack  growth  data  may  have  one  or  more 
fastener  holes  per  specimen.  Some  specimens  may  or  may  not 
be  fatigue  tested  to  failure.  Also,  every  fastener  hole  in 
each  replicate  test  specimen  may  not  contain  a  measurable  fa¬ 
tigue  crack  or  else  the  crack  is  too  small  or  complex  (e.g., 
multiple  crack  origins  and  branching)  for  fractographic  an¬ 
alysis.  A  statistical  scaling  technique  has  been  developed 
[2]  for  determining  the  EXFSD  for  a  "single  hole  population" 
based  on  the  largest  fatigue  crack  per  specimen.  Hence,  it 
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is  necessary  to  read  only  the  fractographic  results  for  the 
largest  crack  per  specimen.  Essential  elements  are  concept¬ 
ually  described  in  Fig.  3-8.  This  technique  is  very  general 
and  is  independent  of  the  distribution  functions  used.  It 
accounts  for  the  number  of  fastener  holes  per  test  specimen 
in  a  given  fractographic  data  sot.  It  minimizes  the  fracto¬ 
graphic  reading  requirements,  permits  a  maximum  utilization 
of  the  available  fractographic  data  and  allows  for  "mixing 
and  matching"  of  fractographic  data  for  the  largest  crack  in 
specimens  with  a  different  number  of  holes. 

Details  of  the  statistical  scaling  technique  developed 
are  given  in  Volume  I  [2].  Essential  features  and  key  aqua¬ 
tions  are  summarized  in  the  following. 


Let  the  cumulative  distribution  of  EIFS  for  a  single 

hole  population  be  denoted  by  Fa(Q)(x)/  and  that  of  the  EIFS 

based  on  the  largest  fatigue  crack  per  specimen  with  i  fastener 

holes  be  denoted  by  F,  ,„.(x).  Assuming  that  fatigue  crack- 

a*  (oj 


ing  in  each  fastener  hole  of  a  specimen  is  statistically 


dependent  of  the  other  holes,  F 
through  the  following. 


a,(0) 


x)  is  related  to  F 


a(0) 


in- 

(x) 


(3-8) 


where  JL  »  number  of  fastener  holes  per  specimen.  Similar  ex¬ 
pressions  for  the  cumulative  distribution  of  TTCI  are  given  in 
Eqs.  (3-9  )and  (3-10j  , 


Vt) 


1 


(3-9) 


FT^(t)  -  1  -  [1  -  FT(t))'t 


(3-10) 


where  PT(t)  »  cumulative  distribution  of  TTCI  for  a  single 
hole  population,  and  F T  {$  ■  cumulative  distribution  of  the 
minimum  TTCI  per  specimen  with  jL  holes. 


In  a  similar  manner,  let  (x)  denote  the  cumulative 
distribution  of  crack  size  at  any  time  t  for  a  single  hole 
population  and  Fa^(t)  M  denote  the  cumulative  distribution 
of  crack  size  at  any  time  t  based  on  the  largest  crack  size 


in  a  specimen  with  JL  holes. 


as  follows. 


Then'  Fa(t)(x) 


is  related  to 


a,  (t) 


(x>  “  CFa(t)Wi 


(3-11) 


The  simple  scaling  technique  described  in  this  section 
has  been  incorporated  into  the  procedure  for  estimating  the 
EIFSD  parameters.  Computer  software  for  the  IBM-compatible 
PC  is  available  for  estimating  the  EIFSD  parameters,  and  for 
checking  goodness-of-fit  [24].  The  validity  of  the  statisti¬ 
cal  scaling  technique  is  evaluated  and  demonstrated  in  Volume 
II  [3]. 


3.3.5  Combined  Least  Square  Sums  Approach  (CLSSA) 

with  the  procedures  described  previously,  we  have  M  EIFS 
data  sets  corresponding  to  M  available  fractographic  data 
sets.  However,  these  M  EIFS  data  sets  are  not  necessarily 
homogeneous,  since  each  data  set  may  have  a  different  scaling 
factor.  For  instance,  the  ith  EIFS  data  set  may  represent 
the  largest  crack  per  specimen  with  ^  holes,  whereas  the  jth 
EIFS  data  set  may  represent  the  largest  crack  per  specimen 
with  /j  holes,  where  In  this  case,  the  ith  and  jth 
EIFS  data  sets  are  non-homogeneous .  Hence,  methods  should  be 
developed  to  utilize  such  data  sets  to  determine  the  EIFSD 
parameters.  In  this  connection,  a  method  has  been  developed 
in  Vol.  I  [2],  referred  to  as  the  "combined  least  square  sums 
approach".  This  approach  enables  one  to  "mix  and  match" 
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these  inhomogeneous  EIFS  data  sets  for  the  determination  of 
EIFSD  parameters.  Step-by-step  procedures  are  described  in 
the  following  for  estimating  the  EIFSD  parameters  of  the  Wei- 
bull-compatible  distribution  (i.e.,  «  and  for  a  given  xu) . 

1.  Define  the  scaling  factor  for  each  EIFS  data  set 
(see  Section  3.3.4).  The  scaling  factor  for  the  ith  EIFS 
data  set  is  denoted  by  j^,  i.e.,  is  the  number  of  fastener 
holes  per  replicate  specimen  in  data  set  i. 

2.  Rank  the  EIFSs  for  each  data  set  separately.  Let 
be  the  jth  smallest  EIFS  value  in  the  ith  data  set. 

3.  Assume  the  Weibull  compatible  distribution  function, 
Eq.  (3-1}  is  used  to  represent  the  EIFSD.  Other  suitable 
distribution  functions  could  be  used,  for  example,  the 
lognormal  compatible,  lognormal  and  two-parameter  Weibull 
distribution  functions. 


4.  Determine  the  Weibull  compatible  EIFSD  parameters  in 
Eq.  (3-1,)  i.e.,  and  <)>  for  an  assumed  xu  that  is  the  EIFS 
upper  bound  limit.  xu  is  chosen  with  the  following 
recommended  constraints.  „  for  clearance-fit  fastener  holes: 
largest  EIFS  value  in  M  data  sets  <  xu  <o.05".  Then,  compute 
<*  and  4>  using  Eqs.  (3-12  )and  0-13),  respectively. 


M  M  i 


,,  *'•  u  lV 

Ml  Ml 


Z»,Z  Z  XJ„-Z  Z^ZZy 


«t  = 


i  *  i  i  =  l  j  - 1 
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In  Eqs.  (3-12  )and  (3-13^/X^,  Y^  are  defined  as  follows# 

Ys 

where#  M  *  number  of  EIFS  data  sets  used  to  determine  the 
EIFSD  parameters,  «  the  total  number  of  EIFS  values  for 
the  ith  data  set  and  j  »  rank  in  ascending  order  of  the  EIFS 
value  in  the  ith  data  set#  i.e.,  j  -  1,  2...,  N^. 

Software  is  available  in  Volume  V  [24],  for  an  IBM  or 
IBM-compatible  PC,  for  implementing  the  CLSSA  described 
above. 

3.3.6  Optimization  of  EIFSD  Parameters  and  Goodness-Of-Fit 

The  Weibull  compatible  EIFSD  parameters  <i.e.#6X,  $6 ,  and 
xu)  in  Eq.  (3-1 )  need  to  be  optimized.  Also#  the  candidate 
EIFSD  should  be  tested  for  "goodness-of-fit".  These  aspects 
are  discussed  in  the  following. 

An  iterative  procedure#  based  on  the  total  standard 
error  (TSE) ,  was  developed  [2]  and  evaluated  [3]  for  opti¬ 
mizing  the  EIFSD  parameters  for  the  Weibull  compatible  dis¬ 
tribution  function.  This  procedure,  conceptually  described 
in  Fig.  3-9,  can  be  used  in  conjunction  with  CLSSA  described 
in  Section  3.3.5  and  can  be  implemented  using  software  file¬ 
name  -  "WCIFQ"  from  Volume  V  [24]. 

The  optimization  procedure  is  as  follows. 

1.  Assume  a  value  for  the  EIFS  upper  bound  limit,  xu# 
within  the  recommended  range  for  clearance-fit  fastener  holes 
(i.e.,  largest  EIFS  value  in  M  data  sets  <  xu  <0.05"). 
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Figure  3-9.  General  Procedure  for  Optimizing  SIFSD  Parameters 
and  Checking  Goodness-of-Fit  for  Compatible  Type 
EIFSD  Func t ion . 
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2.  Determine  <X  and  ^  for  the  assumed  xu  value  using  the 
CLSSA  described  in  Section  3.3.5  and  Eqs.  (3-12 ) and  (3-13) , res¬ 
pectively. 

3.  Compute  the  total  standard  error,  TSE,  for  M  EIFS 
data  sets  using  the  ot ,  £  and  xu  values  from  s  teps  1  and  2 
above  as  follows. 


All  terms  in  Eq.  (3-15) are  the  same  as  those  defined  for  Eqs. 
(3-12  )and  <3-13  J  . 

4.  Repeat  Steps  1-3  for  different  xu  values  and  the 
optimal  (  ot ,  $  ,  xu)  is  obtained  when  the  corresponding  TSE 
is  a  minimum.  Verify  the  goodness-of-fit  for  the  resulting 
EIFSD  using  the  fractographic  data  sets  that  have  been  used  to 
estimate  the  EIFSD  parameters.  For  example,  correlate  theo¬ 
retical  predictions  for  (i)  the  probability  of  crack  exceed¬ 
ance,  p(i,T),  at  a  given  service  time,  T,  and  (ii)  the  cum- 
ultive  distribution  of  service  time  to  reach  any  crack  size 
xl'  FT(x  ) (*) /  with  actual  fractographic  results  for  those 
data  se4s  that  have  been  used  to  define  the  IFQ,  see  Fig.  3~ 
10.  Other  fractographic  data  sets  (e.g.,  for  different 
stress  levels,  load  spectra,  %  bolt  load  transfer,  etc.)  that 
have  not  been  used  to  estimate  the  EIFSD  parameters  can  also 
be  used  to  justify  the  candidate  EIFSD  for  durability  analy¬ 
sis. 


3.4  TWO-SEGMENT  DETERMINISTIC-STOCHASTIC  CRACK  GROWTH 
APPROACH  (DCGA-SCGA) 

The  EIFS  distribution  established  previously  will  serve 
as  a  basis  from  which  the  extent  of  damage  for  a  durability 
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critical  component  at  any  service  time  will  be  pradictad. 
Thia  ia  accomplished  by  growing  tha  antira  EIFS  population 
(or  diatribution)  forward  undar  tha  daaign  loading  apactra. 
In  growing  tha  EIFSD  forward/  tha  computation  of  tha  crack 
growth  damage  accumulation  in  aarvica  ia  dividad  into  two 
aagmenta  for  simplicity.  In  tha  first  aagmant  in  which  the 
crack  aiza  ia  smaller  than  the  rafaranca  crack  size  aQ  for 
TTCI,  a  deterministic  crack  growth  rata  modal/  Eg.  (3-16;,  is 
uaed. 

<*a(t)  -  a(t)  <  aQ 

dt  ;  ' 

In  Eg.  (3-16;  Q1  and  b1  are  empirical-based  constants  de¬ 
pending  on  the  expected  service  loading  spectrum.  They  oan 
ba  obtained  based  on  either  applicable  fractographic  data  or 
crack  growth  predictions  using  a  suitbla  analytical  crack 
growth  program  [e.g.(  31,32].  Since  EZFS  values  are  deter¬ 
mined  from  fractographic  data  in  the  AL-AU  range  by  back-ex¬ 
trapolation  using  the  deterministic  crack  growth  approach,  tha 
EIFS  distribution  must  be  grown  forward  up  to  the  reference 
crack  size  aQ,  that  is  usually  egual  to  AU ,  based  on  tha  de¬ 
terministic  crack  growth  approach. 

The  following  stochasitc  crack  growth  rate  model  is  used 
for  crack  sizes  >  a0 

ckct)/dt  -  X <2 2 [act)]  *.  a.(t)  (3-17) 

in  which  X  is  a  lognormal  random  variable  with  a  median  of 
one;  Q2  and  b2  are  crack  growth  rate  parameters  depending  on 
tha  service  load  spectrum.  Eguat ion (3-17 )  accounts  for  the 
crack  growth  rate  variability  and  is  referred  to  as  the  "log- 
normal  random  variable  modal"  proposed  by  Yang  et  al  [25,26, 
28,33-36] . 

The  probability  density  function  of  the  lognormal  random 
variable  X  with  a  median  1.0  is  given  by 
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(3-18) 


-O 


;  U.2o 

;  u<c 


In  which  <TZ  is  ths  standard  dsviation  of  ths  normal  random 
variable  Z  *  log  X.  Equation  ( 3-18 )  is  used  whan  <rz  is  esti- 
matad  using  the  log  to  base  10  form.  If  i*  based  on  the 
natural  log  form,  fx(u)  given  in  Eq. (3-19) should  be  used. 

(3-19) 

a  O  j  CL  <0 


Note  that  <rz  based  on  the  log  to  base  10  is  equal  to  that 
based  on  the  natural  log  divided  by  the  natural  log  of  10. 
Details  for  estimating  ag  are  given  elsewhere  [2,3]. 


Let  T  be  the  time  for  the  EIFS,  a(G),  to  reach  a  refer¬ 
ence  crack  size  aQ.  Then,  integrating  Eq. (3-16) from  t  -  o  to 
t  «  T  for  b1  *  1,  one  obtains 


r  =  a! £~[at/<u<h 


(3-20) 


in  which  it  is  understood  that  a(T)  -  aQ. 

In  the  region  where  a  (T)  >  aQ  (or  T>  T) ,  Eq.(  3-17)is  in¬ 
tegrated  with  b2  -  1  from  t  -  T  to  t  -  T  (or  from  a(T)  -  aQ 
to  a(t)  -  afT");  with  the  result 


T 3  T'—  J  (3-21) 

Equating  Eqa.  (3-20  )and  0-21)  laada  to  the  following  rela 
tlon  between  a (T )  and  a(0) 


a(o)  =  a.<sxr(-<?lr)l<w!>/aS\  ;<uri>4l 


(3-22) 


In  which 


y  -  (3-23) 

Whan  tha  crack  size,  a (7*),  at  any  sarvlca  time  T  la 
smaller  than  aQ/  tha  ralatlon  batwaan  a (7)  and  a(0)  la  ob- 
tainad  by  intagrating  Eq.  (3-16)  for  bx  -  1  from  t  »  0  to  t  »7* 
aa  follows: 


CUo)  *  4.(T)exr(-Ql7)  ;  <*(l)«Z0  (3-24) 

Equations  (3-22  ) and  (3-24  ) show  tha  ralatlon  batwaan  tha  crack 
size  in  service  a(7“)  and  the  EIFS  a(0).  They  will  be  used  to 
transform  the  EIFSD  to  the  distribution  of  the  crack  size 
a (T)  later. 


3.4.1  Probability  of  Crack  Exceedance  in  Service 

The  probability  that  a  crack  in  the  ith  stress  region 
will  exceed  any  specified  crack  size  xl  at  any  service  time 
7'  is  referred  to  as  the  probability  of  crack  exceedance,  de¬ 
noted  by  p(i,7~).  Depending  on  tha  crack  size  of  interest 
x^,  the  crack  exceedance  probability,  p(i/T),  can  be  derived 
in  the  following  manner. 


3-28 


(1)  When  the  crack  size  of  interest  x1  is  smaller  than 
the  reference  crack  size  aQ,  the  distribution  function 
F  .  (*l)  ■  P[a(T)<x1]  of  the  crack  size,  a (T),  for  x^Sq 
can  be  derived  from  the  distribution  function  of  a(0)  through 
the  transformation  of  (Eq.  3-24; 

fart  - 

in  which 

r)  -  X,  c*?  (-Qt  T) 


(3-25) 


(3-26) 


The  crack  exceedance  probability,  p(i,7),  is  given  by 


fUjri  *  /*[*.(?)  >*,]  -  t  - 

(3-27) 

=  *>&ao 

where  Faf0^(x)  is  the  distribution  function  of  EXFS,  a(0), 
given  by  Eq.  (3-1)  or  other  suitable  distribution  functions. 


(2)  When  the  crack  size  of  interest  x.^  is  larger  than 
aQ,  the  conditional  distribution  function  of  a  (j)  at  any  ser¬ 
vice  time  T  ,  given  X«u,  can  be  derived  from  that  of  a ( 0 ) 
through  the  transformation  of(Eq.  3-22).  Then,  the  uncondi¬ 
tional  distribution  function,  Fft^  (xi)  •  *(“7“)  can  b0  ob~ 
tained  using  the  theorem  of  total  probability;  with  the  re¬ 
sult. 


fart 


(3-28) 


in  which  the  lognormal  probability  density  function  fx(u)  is 
given  by  Eq.  (3-18)  or  (3-19)  and 


<S6f,;rjx-«-)  -  ‘Z()e*r(-4,T)[x,/aJ 


a 


(3-29) 


The  crack  exceedance  probability,  p(i,T),  for  x1>aQ  is 
given  by  p(i,r)  -  1  -  Fa(r)  (x1) ,  i.e.( 

-  /  -/* ft  foj  r|,)( =«.)]•£,«)  ,/u.  (3-30) 

<7 

When  the  Weibull  compatible  distribution,  (Eg.  3-1),  is 
used  for  the  EIFSD,  the  condition  that  Fa(Q) CG(x1jT|x*u) ]  -  1 
for  Gfx.^;  T  |x-u)  >  xy  should  be  reflected  in  the  computer 

program  [24]  for  computing  the  crack  exceedance  probability 
P  ( A 1 7")  1  (Eg  •  3 — 3  0 1) 

3.4.2  Cumulative  Distribution  of  Service  Time 
To  Reach  any  Specified  crack  Size 

Let  TCx^)  be  the  time  for  a  crack  to  reach  any  given 
crack  size  x1  and  FT^X  j (7)  be  the  corresponding  cumulative 
distribution  function,  i.e.,  ET(X  j  (T)  -P[T(x1)<r].  The 
distribution  function  of  T(x^)  is  thil  probability  that  the 
crack  will  reach  a  crack  size  x1  before  service  time  7'. 
Such  a  probability  is  equal  to  the  probability  that  the  crack 
size  a(7“)  at  service  time  will  exceed  x^,  which  is  simply 
the  probability  of  crack  exceedance,  Hence, 

W*  ‘  ^  ^  (3-31) 

=.  rfart  *  *,J  * 

Consequently,  FT ^  (T)  is  obtained  for  any  given  crack  size 
x^  by  computing  the 'crack  exceedance  probability,  p(i,T) ,  at 
different  values  of  service  time  T» 

The  cumulative  distribution  of  service  time,  FT^X  ^  (' T), 
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for  a  crack  to  roach  any  given'  crack  size  x1  ia  determined 
uaing  CEq.  3-314  FT^X  j  (T)  is  obtainad  for  xx  <  aQ  and  for 
xl>a0  by  computing  p?i ,T)  at  diffarant  aarvica  times,  T  ,  ua¬ 
ing  Eq.  (3-27)  and  (3-304  respectively. 

3.5  DURABILITY  ANALYSIS  PROCEDURE 


Tha  durability  analyaia  procadura  for  implementing  tha 
two-aagmant  DCGA-SCGA  daacribad  in  Fig.  3-4  inoludaa  tha  fol¬ 
lowing  baaic  atapa. 

1.  Decide  at  what  level  the  durability  analyaia  will  be 
performed  (e.g.,  single  part,  several  diffarant  parts,  com¬ 
ponent,  etc.). 

2.  Determine  which  structural  details  will  be  included 
in  tha  durability  analyaia  (e.g.,  faatanar  holes,  luge,  cut¬ 
outs,  filiate,  etc.). 

3.  Determine  tha  IFQ  or  suitable  EIFSD  for  each  type  of 
structural  detail  to  be  considered  as  described  in  Section 
3.3. 

4.  For  each  part,  component,  etc.,  group  the  structural 
details  by  type  into  m  stress  regions  where  the  maximum 
stress  in  eaoh  region  may  reasonably  be  assumed  to  be  equal 
for  every  location  or  detail  (e.g.,  fastener  hole). 


5.  Determine  service  crack  growth  parameters  Q1#  Q2  and 
<rz  for  each  stress  region. 

6.  Compute  the  probability  of  crack  exceedance, 
pCi,^),  at  a  given  service  time  T  for  each  stress  region. 


7.  If  desired,  compute  the  cumulative 
service  time  to  reach  a  given  crack  size  x^ 
computation  is  optional. 


distribution  of 

FT<*()«.  This 
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8.  Estimate  the  extant  of  damage  mean  (P  ■  0.5)  and 
upper  bound  limit  for  selected  exceedance  probability  (e.g. , 

?  »  0.05) . 

Essential  equations  and  details  for  implementing  the 
above  steps  are  given  in  the  following  subsections. 

Durability  analysis  guidelines  are  given  in  Section  IV. 
Procedures  and  methods  are  illustrated  in  Section  V  and  else¬ 
where  [2-7,24-29].  Documented  software,  with  a  plotting  cap¬ 
ability,  is  available  in  Volume  V  [24]  for  an  IBM  or  IBM-com¬ 
patible  PC  for  implementing  the  durability  analysis.  This 
software  is  briefly  described  in  Section  VI. 

3.5.1  Service  Crack  Growth  Parameters  Q2  and 

The  relation  between  the  crack  size,  a (f) ,  at  any  ser¬ 
vice  time  T  and  the  EIFS,  a(0) ,  such  as  Eqs.  (3-22  )  and  (3 — 24V 
is  referred  to  as  the  "service  crack  growth  master  curve" 
(SCGMC) .  The  SCGMC  in  each  stress  region  is  determined  by 
either  available  fractographic  results  or  LEFM  crack  growth 
analysis,  in  the  latter  case,  the  LEFM  crack  growth  computer 
program  [e.g.,  31,32]  is  "tuned"  or  "curve-fitted"  to  the 
EIFS  master  curve  in  the  AL-AU  crack  size  region  where  base¬ 
line  fractographic  data  are  available.  Normal  assumptions 
for  the  crack  shape  and  geometry  are  reflected  in  the  crack 
growth  analysis.  Then,  the  SCGMC  is  fitted  by  Eq.  (3-16)  for 
the  crack  size  smaller  than  aQ  to  determine  the  parameters  Q.^ 
and  b1,  and  by  Eq.  (3-17) for  the  crack  size  larger  than  aQ  to 
determine  parameters  Q2  and  b2  (with  X  -  1.0),  using  the 

least  squares  fit  procedure.  The  special  case  b^  ■  b2  -  l.o 
can  be  used  in  Eqs.  (3-16  )and(3-17,)  which  has  been  shown  to  be 
quite  reasonable  for  durability  analyses  [1-3,5-7,17,21,22, 
25,29] . 

If  suitable  fractographic  results  are  available,  it  may 
be  feasible  to  develop  expressions  for  Q1  (Eq.  3-16)  and  Q2 
(Eq.  3-17)  in  terms  of  the  maximum  stress  level  using  an  em- 
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pi'rical  model  proposed  by  Yang  and  Manning  [1-3,5,16]. 


In  Eq.  (3-32,)  -  service  crack  growth  parameter  either  Q1  or 

Q2  for  the  ith  stress  region,  ££  ■  maximum  stress  level  in  the 
ith  stress  region,  and  C  and  V  are  empirical  constants.  Both 
c  and  V  can  be  determined  from  available  base-line  data  or 
suitable  analytical  crack  growth  results  for  Q.  versus  oJ(i  * 
1,  2,  . ..,  m)  using  a  least  square  fit  procedure. 

The  following  notations  are  used  to  determine  the  stan¬ 
dard  deviation,  <j£,  shown  in  Eqs.  (3-18  )and (3-19  )using  a  par¬ 
ticular  fractographic  data  set.  Let  m  *  the  total  number  of 
fatigue  cracks  in  the  fractographic  data  set,  N,,  -  number  of 
da(t)/dts  versus  a(t)s  in  the  AL-AU  range  for  the  jth  fatigue 
crack,  N  «  ■  total  number  of  [da(t)/dt,  a(t)]  pairs  in 

the  AL-AU  range,  (da(t)/dt) ^  "  the  kth  crack  growth  rate 
value  for  the  jth  fatigue  crack,  ajft^)  “  crack  size  for  the 
j th  fatigue  crack  at  the  kth  service  time  t ^  (i.e.,  k  »  1,  2, 
...,  Nj),  Qj  -  crack  growth  rate  parameter  for  the  jth  fa¬ 
tigue  crack  defined  by  Eq.  3-5  and  Q  »  "pooled  Q"  value  for 
the  fractographic  data  set  defined  by  Eq. (3-6  )in  which  Q  - 
* 

The  standard  deviation,  c;,  for  the  ith  stress  region 

z 

shown  in  Eq.  (3-18  )and  (3-19  ) ref  1  acts  the  "log  to  base  10  form" 

and  the  "natural  log  form,"  respectively.  If  the  "log  to 

base  10  form"  is  used,  <r_  is  computed  from  either  Eq.  (3-33) 

z 

or  (3-34.) 

(3-33) 


(3-34) 


JJL. Z  j?  ['h  -^9  A}  a.- UA) ]  4 
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If  the  "natural  log  form"  is  used,  ^  is  computed  using 
either  Eq.  (3-35)  or  (3-36) 

<5  * 


(3-36) 

3.5.2  Probability  of  Crack  Exceedance 

Given  Q1,  Q2  and  cz>  the  probability  of  crack  exceed¬ 
ance,  p(i,  T) ,  for  each  stress  region,  i,  at  any  service 
life, T,  can  be  computed  as  described  in  Section  3.4.1. 

3.5.3  Cumulative  Distribution  of  Service  Time 
To  Reach  Any  Specified  crack  Size 

The  cumulative  distribution  of  service  time  to  reach  any 
crack  size  FT^X  ^(T)#  can  be  computed  for  the  desired 

stress  regions.  This* computation  is  optional  since  the  re¬ 
sults  are  not  needed  for  estimating  the  extent  of  damage. 
The  equations  for  computing  FT^X  j(T)  are  given  in  Section 
3.4.2.  1 

3.5.4  Statistical  Estimation  for  Extent  of  Damage 

The  extent  of  damage  can  be  defined  by  the  statistics  or 
distribution  of  the  number  of  structural  details  or  ligaments 
in  the  durability-critical  component  expected  to  exceed  spe¬ 
cified  crack  size  limits  at  a  given  service  time.  From  a 
functional  impairment  standpoint,  the  extent  of  damage  may  be 
interpreted  as  the  number  of  locations  where  the  accumulated 
crack  size  exceeds  limiting  crack  sizes  for  functional  im¬ 
pairment.  For  example,  a  through- the-thickness  crack  in  a 
fuel  tank  nay  cause  fuel  leakage  and  the  dimension  between 
adjacent  structural  details  may  be  considered  as  a  crack  size 
limit  for  ligament  breakage.  The  mean  and  upper  bound  limit 


(3-35) 
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(see  Fig.  3-5)  for  the  extent  of  damage  can  be  estimated  for 
selected  exceedance  probabilities  as  follows. 

The  number  of  details,  N(>i.  ,T  ),  in  the  ith  stress  re¬ 
gion  with  a  crack  size  greater  than  at  the  service  time  T, 
is  a  statistical  variable.  The  mean  value,  N (X,T),  and  the 
standard  deviation,  <3^ (i,  T )  ,  are  determined  using  the  Bino¬ 
mial  distribution  [50]. 

J7U/r)  ^  a*  pUj  r)  (3-37) 

(3-38) 

in  which  denotes  the  total  number  of  details  in  the  ith 
stress  region.  The  average  number  of  details  with  a  crack 
size  exceeding  x1  at  the  service  time  T  for  m  stress  regions, 
C(T) ,  and  the  standard  deviation,  <£(r) ,  can  be  computed  us¬ 
ing  Eqs.  (3-39)  and (3-40;  respectively. 


Lrt  *  £/7?4r) 

A  —! 

(3-39) 

rr.-n  - 

(3-40) 

Equations  0-39)  and  (3-40) can  be  used  to  quantify  the  extent  of 
damage  for  a  single  detail,  a  group  of  details,  a  part,  a 
component,  or  an  airframe.  L (T)  approximately  corresponds 
to  a  50%  probability.  Upper  and  lower  bound  limits  for  the 
"extent  of  damage"  can  be  estimated  using  the  Binomial  dis¬ 
tribution,  e.g.,  L (T)  +  Z(J^(T) ,  with  Z  being  the  number  of 
standard  deviations,  from  the  mean,  L(T ) .  For  example,  Z  ■ 
1.65  and  Z  »  -1.65  correspond  to  exceedance  probabilities  of 
P  -  0.05  and  0.95,  respectively.  Equations  (3-37  )to  (3-40)  are 
valid  assuming  that  the  crack  growth  accumulation  for  each 
detail  is  statistically  independent  [50]. 
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SECTION  IV 


DURABILITY  ANALYSIS  GUIDELINES 

Durability  analysis  guidalines  are  presented  in  this 

section  for  the  following:  (1)  acquisition  and  utilization  of 

fractographic  data,  (2)  determination  of  suitable  EIFSD  for 

durability  analysis,  (-'  determination  of  service  crack 

growth  master  curves  fv  ^all  and  large  crack  size  regions, 

(4)  estimation  of  standard  deviation,  CL  for  the  large  crack 

z 

size  region,  (5)  initial  flaw  size  considerations,  and  (6) 
extent  of  damage.  Practical  aspects  are  emphasized. 


4.1  ACQUISITION  AND  UTILIZATION  OF  FRACTOGRAPHIC  DATA 

Guidelines  are  presented  in  this  section  for:  (1)  test 
specimens,  (2)  fatigue  testing,  (3)  fractography,  (4)  utiliz¬ 
ing  existing  fractographic  data  and  (5)  screening/censoring 
fractographic  data. 

Initial  fatigue  quality  (IFQ)  data  can  be  acquired  for 
various  materials  and  structural  details  as  a  part  of  the 
Aircraft  Structural  Integrity  Program  (ASIP)  effort.  Speci¬ 
mens  tested  under  ASIP  can  provide  data  applicable  to  both 
"durability"  and  "damage  tolerance".  For  example,  if  struc¬ 
tural  details  in  test  specimens  are  not  preflawed,  "natural 
fatigue  crack"  data  can  be  used  not  only  to  estimate  the  IFQ 
of  structural  details  but  also  to  satisfy  durability  and  dam¬ 
age  tolerance  data  requirements.  Depending  on  the  degree  of 
confidence  desired  and  circumstances,  additional  tests  and 
fractographic  evaluations  may  be  desirable  to  estimate  the 
IFQ.  In  the  future,  ASIP  test  plans  should  be  designed  to 
satisfy  the  needs  for  initial  fatigue  quality,  durability  and 
damage  tolerance.  This  approach  will  minimize  the  IFQ  data 
acquisition  costs  with  a  minimal  impact  on  schedule. 
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4.1.1  Test  Specimens 


Initial  fatigue  quality  (IFQ)  data  can  be  acquired  from 
suitable  fractographic  results  for  the  type  of  structural  de¬ 
tail  to  be  reflected  in  the  durability  analysis  (e.g.,  fast¬ 
ener  holes,  cutouts,  lugs,  fillets,  etc.).  So  far,  IFQ  of 
clearance-fit  straight-bore  and  countersunk  fastener  holes 
has  been  primarily  investigated  [3,5-7,17,30-22,27,28,30]. 
Such  investigations  were  concerned  with  clearance-fit  fasten¬ 
ers  in  holes  without  special  life  enhancement  features  (e.g., 
cold  working,  interference  fit  fasteners,  force-fit  bushings, 
etc . ) . 


The  test  specimen (s)  used  to  acquire  the  IFQ  data  should 
include  the  type  of  structural  detail  for  which  the  IFQ  is 
sought.  Also,  the  structural  detail  should  reflect  the  ap¬ 
plicable  manufacturing,  assembly  and  processing  methods,  in¬ 
cluding  fastener  hole  life  enhancement  considerations.  Fur¬ 
ther  research  is  needed  to  develop/evaluate  appropriate  test 
specimen  designs  for  acquiring  IFQ  data  for  cutouts,  fillets, 
lugs,  etc. 

Various  types  of  specimen  have  been  previously  used  to 
acquire  IFQ  data  for  clearance-fit  fastener  holes  [e.g.,  37- 
39].  Three  such  specimen  configurations  are  shown  in  Fig. 
4-1.  Here  are  a  few  comments  about  these  specimen  based  on 
our  experience.  The  simple  dog-bone  specimen  shown  in  Fig. 
4-1 (a)  with  a  single  fastener  hole  and  a  bolt  installed  is 
best  suited  for  acquiring  the  IFQ  data.  The  double-reversed 
dog-bone  specimen  in  Fig.  4-1 (b)  is  simple  for  fatigue  test 
to  acquire  the  fractographic  data.  However,  there  is  no  way 
to  accurately  control  the  amount  of  bolt  load  transfer  since 
the  percent  of  bolt  load  transfer  depends  on  the  fastener- 
hole  fit  and  the  applied  load  level.  A  multi-hole  type  speci¬ 
men,  such  as  that  shown  in  Fig.  4-l(c),  is  potentially 
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(a)  No  Bolt  Load  Transfer  Specimen 


(b)  Double  Reversed  Dog-Bone  Specimen  Design 


(c)  Multi-Hole  No  Bolt  Load  Transfer  Design 


Figure  4-1.  Common  Types  of  Specimens  Used  for  Acquirin 
IFQ  Data  for  Fastener  Holes. 


attractive  for  acquiring  the  IFQ  data,  however,  it  has  the 
following  limitations  or  shortcomings. 

Fatigue  cracks  may  initiate  sooner  in  one  hole  than  an¬ 
other.  Ideally,  one  may  assume  that  each  fastener  hole  in 
the  specimen  of  Fig.  4-1 (c)  is  equally  stressed  until  failure 
occurs.  However,  our  experimental  results  [37-39]  show  that 
each  hole  may  not  be  equally  stressed  during  spectrum  fatigue 
testing  although  the  fastener  holes  are  spaced  far  enough 
apart  to  minimize  the  effects  of  adjacent  holes  on  the  crack 
initiation  and  crack  growth  in  neighboring  holes.  There  is 
apparently  some  interaction  between  the  fatigue  cracks  in 
adjacent  holes  for  the  larger  fatigue  cracks.  The  three-hole 
specimen  designs  shown  in  Fig.  4-1 (c)  may  be  reasonable  for 
acquiring  small  fatigue  crack  data,  but  it  is  not  recommended 
for  acquiring  large  fatigue  crack  growth  results  in  which 
there  may  be  some  interaction  effect  between  adjacent  holes. 

Three  types  of  specimen  are  recommended  for  acquiring 
IFQ  data  for  fastener  holes  as  shown  in  Fig.  4-2.  In  Fig. 
4-2 (a)  the  simple  dog-bone  specimen  with  a  single  fastener 
hole  is  well  suited  for  acquiring  IFQ  data.  The  two-for-one 
specimen  shown  in  Fig.  4-2 (b)  is  also  attractive  for  acquir¬ 
ing  fatigue  cracking  data.  For  example,  the  test  specimen 
can  be  fatigue  tested  to  failure  through  one  fastener  hole. 
Then,  the  specimen  is  reworked  by  sawing  off  the  broken  piece 
and  making  an  end  lug  cut  of  the  center  built  up  area.  Test¬ 
ing  is  continued  until  failure  occurs  in  the  second  hole. 
When  bolt  /load  transfer  is  to  be  accounted  for  in  the  IFQ 
representation  of  the  specimen,  the  specimen  design  shown  in 
Fig.  4-2 (c)  is  well  suited.  With  this  design,  the  basic  spec¬ 
imen  can  be  tha  same  as  that  shown  in  Fig.  4-2 (a).  The  only 
difference  is  that  a  loading  bar  is  used  to  directly  load  the 
bolt  in  the  center  hole  (either  in  double  or  single  shear) . 
Both  the  lug  end  and  the  loading  bar  are  connected  to  sepa¬ 
rate  loading  rams  which  are  synchronized  for  spectrum  fatigue 
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(a)  No  Bolt  Load  Transfer  Dog-Bone  Specimen 


<b) 


Two-Far-One  No-Bolt  Load  Transfer  Dog-Bone 
Specimen 


(Use  double  or  single  shear  configuration) 


(c)  Dog-Bone  Specimen  With  Bolt  Load  Transfer 


Figure  4-2. 


Recommended  Specimen  T^ 
IFQ  Data. 


>es  for  Acquiring 
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testing. 


The  biggest  shortcoming  of  the  double-reversed  dog-bone 
specimen  (see  Fig.  4-1 (b) )  for  acquiring  IFQ  data  under  bolt 
load  transfer  conditions  is  that  the  amount  of  bolt  load 
transfer  varies  and  cannot  be  controlled.  The  recommended 
bolt  load  transfer  specimen  design  shown  in  Fig.  4-2 (c)  re¬ 
quires  a  more  complicated  test  setup  and  more  expensive  test¬ 
ing  costs.  However,  the  dog-bone  specimen  shown  in  Fig.  4- 
2(a)  can  be  used  for  the  setup  shown  in  Fig.  4-2 (c).  What¬ 
ever  specimen  design  is  used  to  acquire  IFQ  data  under  bolt 
load  transfer  conditions,  the  amount  of  bolt  load  transfer 
should  be  controlled  so  that  the  effect  of  load  transfer  on 
IFQ  can  be  more  readily  defined. 

Comments  on  test  specimens  geometry  and  recommended  num¬ 
ber  of  test  specimens  for  acquiring  the  IFQ  data  are  as  fol¬ 
lows.  If  possible,  the  test  specimens  should  be  vide  enough 
to  acquire  valid  crack  growth  data  up  to  a  1"  crack  size. 
The  IFQ  is  estimated  using  fractographic  data  in  the  small 
crack  size  region  (e.g,  AL-AU  -  0.01"  -  0.05").  However,  the 
large  crack  size  data  is  very  useful  for:  (1)  justifying  the 
candidate  EIFSD  for  applications  in  large  crack  size  region, 
and  (2)  developing  suitable  service  crack  growth  master 
curves  for  desired  durability  analysis  conditions.  The  spec¬ 
imen  cross  section  and  pin-to-pin  length  should  provide  a 
stable  specimen  under  applicable  compressive  loads  in  the 
spectrum. 

If  wide  test  specimens  are  not  practical  due  to  severe 
time  and  cost  constraints,  narrow  width  specimens  will  be  the 
second  choice,  in  which  case  specimens  should  be  tested  to 
failure  under  spectrum  loading.  Fractographic  data  for  the 
small  crack  size  region  (e.g.,  AL-AU  ■  .01"  -  .05")  can  be 
efficiently  acquired  in  this  manner.  This  information  can  be 
used  to  estimate  the  EIFSD  parameters  and  to  justify  the  can- 


4-6 


didate  EIFSD  for  durability  analysis.  However,  the  biggest 
drawback  is  that  the  crack  growth  data  in  the  large  crack 
size  region  are  not  available. 

The  number  of  test  specimens  used  to  acquire  the  IFQ 
data  depends  on  the  cost  constraint,  desirable  confidence 
level  and  time  schedule.  Of  course,  the  more  specimens  the 
better  confidence  in  the  results.  However,  as  a  general 
guideline  the  following  is  recommended.  Use  10-30  specimens 
to  acquire  the  IFQ  data  for  each  test  condition  to  be  consid¬ 
ered. 


4.1.2  Testing  Guidelines 

All  fatigue  tests  for  IFQ  specimens  are  recommended  to 
be  performed  at  room  temperature  in  a  lab  air  environment  un¬ 
til  failure  occurs.  A  test  system  which  automatically  shuts 
off  when  the  specimen  in  the  fixture  fails  is  recommended. 
This  way,  fatigue  testing  can  be  periodically  monitored  and 
new  specimens  can  be  loaded  for  subsequent  testing. 

If  possible,  IFQ  data  should  be  acquired  for  replicate 
specimens  at  three  different  stress  levels  for  a  given  load 
spectrum.  A  high,  low  and  intermediate  stress  level  are 
recommended  to  acquire  the  fatigue  cracking  data  base  where 
possible.  This  type  of  information  is  very  useful  for  deter¬ 
mining  the  EIFSD  parameters,  and  it  provides  a  basis  for 
tuning  the  analytical  crack  growth  program  to  make  durability 
analysis  predictions  for  different  stress  levels.  Further¬ 
more,  this  information  can  be  used  to  establish  an  empirical 
relationship  for  crack  growth  parameter  as  a  function  of 
stress  level, (Eg.  3-32}. 

We  recommend  that  all  specimens  be  tested  to  failure  to 
assure  as  much  uniformity  as  possible  in  the  crack  growth  re- 
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suits  (i.e.,  acquire  fractographic  data  that  cover  a  specifi¬ 
ed  AL-AU  range) .  If  fatigue  testing  is  stopped  at  a  specifi¬ 
ed  tine  rather  than  at  specimen  failure,  the  resulting  fa¬ 
tigue  cracks  in  the  fastener  holes  may  be  too  small  to  reli¬ 
ably  read  the  fractographic  data.  Likewise,  if  the  specimen 
is  not  tested  to  failure,  the  resulting  fatigue  crack,  if 
present,  may  not  be  obvious,  thus,  the  specimen  must  be  bro¬ 
ken  open  to  reveal  the  fracture  surface (s) .  This  requires 
careful  specimen  preparation  to  avoid  possible  damage  to  the 
fracture  surface.  When  the  specimen  is  tested  to  failure,  on 
the  other  hand,  a  "clean"  fracture  surface  results. 


4.1.3  Fractographic  Data  Guidelines 

The  minimum  fractographic  crack  size  to  be  read  should, 
if  possible,  correspond  to  the  crack  size  lower  limit,  AL,  to 
be  used  to  define  the  IFQ.  For  example,  AL  *  0.01"  is  recom¬ 
mended.  If  possible,  the  fractography  should  be  read  to 
cover  the  selected  AL-AU  range  for  defining  the  IFQ  (e.g., 
AL-AU  -  .01"  -  .05").  When  the  fractographic  data  acquired 
does  not  cover  the  selected  AL-AU  range,  TTCIs  for  a  selected 
reference  crack  crack  size,  aQ,  may  have  to  be  extrapolated. 
"Mixing  and  matching"  TTCIs  based  on  extrapolations  and  in¬ 
terpolations  should  be  avoided  where  possible,  because  extra¬ 
polated  values  cannot  be  verified  and  neither  can  the  contri¬ 
bution  to  the  total  variance  of  the  data  base. 

Use  automated  crack  monitoring  techniques  as  much  as 
possible  to  minimize  fractographic  acquisition  costs.  Also, 
automatic  storing  of  the  fractographic  results  directly  into 
the  computer  can  minimize  the  time  and  costs  for  plotting  re¬ 
sults  and  for  estimating  the  EIFSD  parameters. 
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4.1.4  Utilization  of  Existing  Fractographic  Data 


Existing  fractographic  data  should  be  utilized  where 
possible  to  estimate  the  IFQ.  In  some  cases,  test  specimens 
may  contain  multiple  fastener  holes  but  fractographic  results 
may  be  available  for  only  the  largest  fatigue  crack  per  spec¬ 
imen.  In  other  cases,  test  specimen  may  have  only  a  single 
fastener  hole.  Furthermore,  it  may  be  necessary  to  utilize 
fractographic  data  from  different  sources  to  estimate  the 
IFQ.  In  this  connection,  the  statistical  scaling  procedure 
described  in  Section  3.3.4  and  elsewhere  [2,3]  can  be  used  to 
normalize  the  fractographic  results  to  a  single  hole  basis. 


Fractographic  data  used  to  define  the  IFQ  should  first 
be  screened  and  censored  for  data  anomalies  such  as  crack 
growth  data  extremes  (e.g.,  very  fast  or  very  slow  compared 
to  most  results  in  data  sets) .  Software  is  available  in  Vol¬ 
ume  V  [24]  for  efficiently  plotting  and  displaying  the  frac¬ 
tographic  data  for  desired  crack  size  ranges.  This  software 
can  be  implemented  on  an  IBM  or  IBM-compatible  personal  com¬ 
puter.  Screening  and  censoring  is  essential  to  assure  that 
the  data  is  reasonably  homogeneous  and  it  covers  the  selected 
AL-AU  range  to  be  used  to  define  the  IFQ.  "Data  sparsity" 
results  when  a  given  fatigue  crack  has  little  or  no  data  in 
the  desired  AL-AU  crack  size  range. 

The  fractographic  data  base  should  reflect  only  those 
fatigue  cracks  which  originate  in  the  bore  of  the  fastener 
hole.  Ail  surface  fatigue  cracks  should  be  excluded  because 
these  are  not  typical  of  the  fatigue  cracking  process  to  be 
modeled. 


4.2  DETERMINATION  OF  SUITABLE  EIFSD 
The  following  guidelines  apply  to  the  determination  of  a 
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suitable  EIFSD  for  durability  analysis. 

1.  The  EIFS  upper  bound  limit,  xy,  has  the  following 
constraint:  largest  EIFS  in  any  data  set  <  xu  <  0.05".  In 
general,  an  xu  value  between  0.03"  and  0.05"  is  recommended. 
The  maximum  xu  of  0.05"  is  set  by  NDI  considerations,  current 
damage  tolerance  initial  flaw  size  requirements  and  the  eco¬ 
nomical  repair  limit  for  fastener  holes.  If  xu  >  0.05"  is 
allowed,  this  means  that  the  probability  of  exceeding  the 
economical  repair  limit  will  be  greater  than  zero  at  time 
zero.  This  may  not  be  a  realistic  condition  for  newly  manu¬ 
factured  fastener  holes. 

2.  A  fractographic  crack  size  range,  AL-AU  ■  0.01"  - 
0.05"  is  recommended  for  determining  the  IFQ.  Other  cracic 
size  ranges  could  also  be  used.  Whatever  AL-AU  limits  are 
used,  the  pooled  Q  value  for  each  data  set  should  be  deter¬ 
mined  from  the  same  AL-AU  range.  Ideally,  the  fractographic 
data  should  cover  the  selected  AL-AU  range. 

3.  For  AL-AU  *  0.01"  -  0.05",  a  crack  size  of  0.05"  for 
TTCls  and  a  reference  crack  size  of  aQ  -  0.05"  are  considered 
reasonable.  In  any  case,  we  recommend  the  following  limits 
for  defining  aQ:  AL  <  aQ  <  AU. 

4.  We  recommend  that  the  parameters  for  the  Weibull- 
compatible  distribution  function  (i.e.,0d,£  and  xu)  be  esti¬ 
mated  using  the  "EIFS  fit",  the  combined  least  square  sums 
approach  (CLSSA) ,  data  pooling  and  statistical  scaling  proce¬ 
dure  [2,3].  It  has  been  shown  that  the  same  EIFSD  parameter 
values  can  be  obtained  using  either  the  "EIFS  fit"  or  the 
"TTCI  fit"  when  the  CLSSA  is  used  [3]. 

5.  While  values  of  ed  and  $  depend  upon  xu,  we  have 
found,  in  general,  that  comparable  predictions  for  p(i,  T  ) 
and/or  f  (t)  can  be  obtained  using  different  xu  values 

T (X.  ; 
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and  the  corresponding  and  <f>  .  This  is  due  to  the  fact  that*, 

and  $  are  obtained  for  a  given  xu  using  the  least  square  fit 
procedure . 


6.  The  Weibull  compatible  EIFSD  function  given  in  Eq. 
(3-1)  is  recommended  for  defining  IFQ.  Other  distribution 
functions  could  also  be  used  for  this  purpose,  in  any  case, 
a  compatible  type  EIFSD  is  recommended  for  defining  IFQ  be¬ 
cause  such  a  distribution  imposes  physically  meaningful  lim¬ 
its  on  TTCIs  and  EIFSs.  For  example,  all  TTCIs  are  non-neg¬ 
ative  quantities  with  values  >  0  and  the  maximum  EIFS  allowed 
is  governed  by  xu.  The  lognormal  compatible  is  considered  to 
be  another  reasonable  EIFSD  function. 


7.  The  methods  developed  for  determining  the  IFQ  of 
fastener  holes  has  been  evaluated  for  both  straight-bore  and 
countersunk  clearance-fit  fastener  holes.  Further  research 
is  needed  to  account  for  the  effects  of  hole  life-enhancement 
features,  such  as  cold  working,  interference  fit  fastener, 
press-fit  bushings,  etc.  on  IFQ. 


8.  The  candidate  EIFSD  should  be  justified  for  the  plan¬ 
ned  durability  analysis.  As  a  minimum,  the  candidate  EIFSD 
should  be  grown  forward  using  the  applicable  EIFS  master 
curve  for  each  data  set  to  predict:  (1)  p(i,T)  at  a  selected 
service  time,  and  (2)  ft(x  ) ^  at  x1  -  aQ.  Predictions 
should  be  correlated  with  the  ranked  results  for  crack  size 
and  TTCI ,  respectively.  Software  is  available  in  Volume  V 
[24]  for  testing  the  candidate  EIFSD.  Other  fractographic 
data  sets  not  included  in  the  determination  of  the  EIFSD  pa¬ 
rameters  could  also  be  used  to  test  the  candidate  EIFSD  in 
the  small  and  large  crack  size  regions.  If  reasonable  corre¬ 
lations  are  obtained  for  the  candidate  EIFSD  in  the  manner 
described  above,  the  EIFSD  may  be  justified  for  further  dur¬ 
ability  analysis  (e.g.,  same  type  of  load  spectrum  (fighter, 
bomber  or  transport) ,  different  stress  levels  and  %  bolt  load 
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transfer) .  An  EIFSD  is  unacceptable  if  reasonable  correla¬ 
tion  cannot  be  obtained  in  the  areas  of  most  interest,  when 
the  EIFSD  is  grown  forward. 

4.3  ESTIMATION  OF  STANDARD  DEVIATION  C£ 

The  standard  deviation,  <f  ,  of  the  log  crack  growth 
rate  in  the  large  crack  size  region  is  needed  to  implement 
the  two-segment  DCGA-SCGA  described  in  Section  3.4.  If  suit¬ 
able  fractographic  results  are  available,  a  can  be  determin- 

z 

ed  using  Eq.  (3-36  )as  described  in  Section  3.5.2,  otherwise  <jl 
will  have  to  be  assumed.  Ranges  of  °z  value  are  shown  in 
Table  4-1  for  both  countersunk  and  straight-bore  fastener 
holes.  These  results  are  based  on  extensive  fractographic 
data  evaluations  for  7474-T7351  aluminum  with  clearance-fit 
fastener  holes.  The  results  in  Table  4-1  reflect  the  natural 
log  base  and  provide  information  for  a  reasonable  assumption 
for  <*2  value. 


4.4  SERVICE  CRACK  GROWTH  MASTER  CURVE 

A  service  crack  growth  master  curve  (SCGMC)  for  each 
stress  region  is  needed  to  grow  the  EIFSD  forward  to  predict 
p(i,  T  )  at  a  given  service  time,  T,  orp_,v  .  (t)  for  a  given 
crack  size  x^.  Recommended  procedures  and  guidelines  are 
presented  in  the  following  for  determining  a  suitable  SCGMC 
for  the  small  and  large  crack  size  regions. 

4.4.1  Small  Crack  Size  Region 

A  SCGMC  is  needed  for  the  small  crack  size  region  (e.g. 
a(t)  <  0.05”)  to  evaluate  functional  impairment  due  to  exces¬ 
sive  cracking.  In  most  durability  design  situations,  a  suit¬ 
able  LEFM  analytical  crack  growth  program  [e.g.,  31,32]  is 
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Table  4-1  a"z  Range*  for  7475-T7351  Aluainua  for 

Straight-Bore  and  Countersunk  Fastener  Holes. 


Type 

Hole 

range(3) 

m 

1— 

Notes: 


(1)  Straight-bore  fastener  hole  (clearance-fit) 

(2)  Countersunk  fastener  hole  (clearance-fit) 

(3)  AL-AU  *  0 . 05"  -  2" 

(4)  Ref.  Volume  II  (3] 
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used  to  develop  the  SCGMC  for  the  desired  analysis  conditions 
because  applicable  fractographic  results  may  not  be  avail¬ 
able.  The  following  general  procedure  for  developing  a  SCGMC 
is  recommended  for  durability  analysis  applications  in  the 
small  crack  size  region. 

1.  Define  the  bases  for  the  EIFSD  to  be  used  in  the 
durc^ility  analysis  and  for  EIFS  master  curve.  For  example: 
(1)  what  fractographic  crack  size  range,  AL-AU,  was  used?; 
and  (2)  what  method  was  used  to  define  the  EIFS  master  curve, 
including  criterion  for  goodness-of-fit  and  crack  shape? 

Note:  An  empirical  EIFS-Service  time  relationship 
(e.g.,  Eq.  (3-4)) is  recommended  for  general 
applications  so  that  consistent  EIFSs  will 
be  obtained  by  different  aerospace  contrac¬ 
tors  for  the  same  fractographic  data  base. 

2.  use  <?.  suitable  analytical  crack  growth  program  to 
"curve  fir.'*  or  "tune  to"  the  EIFS  master  curve  or  curves  in 
t'-n  £-«?<.. togruphic  crack  size  range,  AL-AU.  The  "curve  fit"  to 
til";  LjCFS  master  curve  in  the  selected  AL-AU  range  is  accom¬ 
plished  using  the  applicable  conditions  reflected  in  the  EIFS 
master  rurve  (i .e.,  load  spectrum,  stress  level,  %  bolt  load 
transfer,  hole  type/diameter) „  The  procedures  are  given  as 
follows:  (1)  pj.ot  the  EIFS  master  curve  to  cover  the  appli¬ 
cable  AL-AU  range;  (2j  select  a  crack  growth  model  (e.g., 
Walker-  uK,  Forman,  r. i:c.);  (3)  select  da/dN  versus  AK  data 
and  calibrate  the  crack  growth  model  parameters  for  given  ma¬ 
larial;  (4)  select  a  loud-retardation  model  (e.g.,  modified 
Wil.lenborg,  Wheeler,  etc.);  and  (5)  by  trial  and  error,  de- 

rmino  the  remaining  model  parameters  required  to  obtain  a 
u;jufu.‘.e  c-. , -■*  fit  (to  the  EIFS  master  curve  in  the  AL-AU 
■ti.r.ge)  using  the  analytical  crack  growth  ■  .  igram.  The  good- 
..»ct-of-f ft  is  determined  subjectivsly.  Ti  Am.  above  procedure 
if  conceptually  described  in  Fig.  4-3 (a). 
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3.  The  tuned  analytical  crack  growth  program  is  then 
used  to  predict  the  crack  growth  over  the  applicable  AL-AU 
range  using  the  applicable  durability  analysis  conditions. 
For  example,  a  specific  spectrum,  stress  level  and  %  bolt 
load  transfer  and  assumed  crack  shape/geometry  are  used  to 
predict  the  crack  size  at  a  given  time.  This  step  is  con¬ 
ceptually  illustrated  in  Fig.  4-3 (b)  for  three  different 
stress  levels  (  c^  >  <r2  >  The  crack  size-time  pre¬ 
dictions  in  the  AL-AU  range  are  indicated  in  Fig.  4-3.  The 
analytical  crack  growth  program  is  used  further  to  make  crack 
size-time  predictions  for  crack  size  greater  than  AU  where 
LEFM  principles  apply.  Procedures  and  assumptions  for  the 
analytical  crack  growth  analysis  are  the  same  as  those  used 
for  a  typical  damage  tolerance  analysis. 

4.  Estimate  Q1  in  Eg.  (3-16  lusing  the  predicted  crack 
size-time  predictions  (i.e.,  a{t)  versus  t)  in  the  AL-AU 
range  for  a  given  stress  region,  depicted  in  Fig.  4-3 (b). 
Methods  for  estimating  are  given  in  Section  III.  As 
shown  in  Fig.  4-3 (b) ,  LEFM  principles  are  used  only  for  the 
crack  size  range  where  such  principles  apply. 


4.4.2  Large  Crack  Size  Region 

For  the  two-segment  DCGA-SCGA,  the  first  segment  covers 
the  small  creek  size  range  (i.e.,  a(t)  <  aQ)  and  the  second 
segment  covers  the  large  crack  size  range,  a(t)  >  aQ,  as 
shown  in  Fig.  4-4.  The  first  segment  was  obtained  previously 
whereas  the  second  segment  can  be  determined  as  follows. 

The  tuned  analytical  crack  growth  program  obtained  above 
is  used  to  predict  tha  crack  growth  (i.e.,  a(t)  versus  t) 
from  aQ  -  AU  to  AU#  for  the  desired  analysis  conditions 
(i.e.,  stress  level,  load  spectrum,  %  bolt  load  transfer, 
etc.).  For  example,  such  crack  growth  predictions  are  shown 
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Figure  4-4.  Uae  of  Analytical  Crack  Growth  Predict Iona  and 
Crack  Growth  Model  to  Establish  SCGhCa  for 
Durability  Analysis  Extension. 


in  Fig.  4-4  for  three  different  stress  levels,  in  which  0^ 

>  <T2  >  In  Fig.  4-4  the  analytical  crack  growth  pre¬ 

dictions  cover  a  crack  size  range  where  LEFM  principles  apply 

(i.e.,  a(t)  >  dg  *  AU) .  The  crack  growth  rate  parameter  Q2 
for  the  second  segment  of  SCGMC  in  Eq. (3-17)  can  be  estimated 
similarly. 

The  two  crack  growth  master  curvj  segments  for  a  given 
stress  level  can  be  physically  combined  into  a  single  SCGMC 
as  illustrated  in  Fig.  4-5.  At  point  i  segments  1  and  2 
have  the  same  (a(t),  t)  values  but  not  necessarily  the  same 
slopes. 

The  two-segment  SCGMC  for  the  DCGA-DCGA  can  also  be  used 
for  the  DCGA-SCGA  for  the  same  crack  growth  analysis  condi¬ 
tions  (e.g.,  stress  level,  load  spectra,  etc.).  The  only 
difference  is  that  the  parameter  <rz  is  required  to  implement 
the  DCGA-SCGA.  is  the  standard  deviation  of  da/dt  with 

respect  to  the  plot  of  In  da/dt  -  In  a(t)  +  Qt.  <TZ  can  be 
determined  from  suitable  fractographic  results,  if  available. 

4.5  GUIDELINES  FOR  SELECTING  STRESS  REGIONS 
AND  STRESS  LEVELS 

For  durability  analysis  purposes,  a  durbility-critical 
part  or  component  is  divided  into  many  stress  regions.  In 
each  stress  region,  the  stress  level  at  each  structural  de¬ 
tail  is  approximately  the  same.  The  number  of  stress  regions 
needed  for  the  durability-critical  component  depends  on  the 
types  of  structural  details  to  be  considered  (e.g.,  fastener 
holes,  cutouts,  lugs,  fillets,  etc.)  and  the  variation  of 
governing  stress  levels.  Different  types  of  structural  de¬ 
tails  cannot  be  included  in  the  same  stress  region  and  they 
should  be  separated.  Some  structural  details  introduce  a  high¬ 
er  stress  intensity  factor  and  hence  a  higher  crack  growth 
rate. 
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Likewise,  appropriate  finite  element  grid  sizes  should  be 
used  to  achieve  the  stress  analysis  accuracy  desired.  A 
suitable  stress  analysis  is  very  important  because  the  gover¬ 
ning  stress  for  a  given  stress  region  can  have  a  significant 
influence  on  the  crack  growth  predictions  for  the  structural 
details. 


4.6  EXTENT  OF  DAMAGE  GUIDELINES 

The  durability  analysis  methodology  developed  can  be 
used  to  estimate  the  extent  of  damage  in  a  durability-criti¬ 
cal  component  due  to  excessive  cracking,  fuel  leaks  and  liga¬ 
ment  breakage.  The  extent  of  damage  depends  on  the  specified 
crack  size  limits  for  functional  impairment.  For  example, 
typical  limits  are  illustrated  in  Table  4-2  for  fastener 
holes.  Functional  impairment  crack  size  limits  for  other 
types  of  structural  details  are  specified  by  the  user. 

Structural  details  in  the  durability-critical  component 
to  be  reflected  in  the  durability  analysis  are  divided  into 
stress  regions.  Functional  impairment  crack  size  limits  can 
vary  for  each  stress  region  and  each  type  of  structural  de¬ 
tail.  The  extent  of  damage  for  a  given  exceedance  probabil¬ 
ity  can  be  predicted  separately  for  each  type  of  structural 
detail  in  the  component  (e.g.,  fastener  holes,  cutouts,  fil¬ 
lets,  lugs,  etc.)  and  each  type  of  functional  impairment. 
Also,  the  overall  extent  of  damage  due  to  different  types  of 
structural  details  can  be  estimated  by  combining  the  extent 
of  damage  results  for  all  details  in  each  stress  region. 

The  extent  of  damage  at  a  given  service  time  for  a  dur¬ 
ability-critical  component  should  be  estimated  for  selected 
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exceedance  probabilities.  For  example,  when  the  binomial 
distribution  is  approximated  by  the  Normal  distribution,  the 
average  extent  of  damage  corresponds  to  an  exceedance  probab¬ 
ility  of  P  -  0.5.  The  upper  bound  limit  for  the  extent  of 
damage  could  be  estimated,  for  example,  at  P  ■  0.05.  In 
other  words,  the  probability  of  exceeding  the  upper  bound 
limit  for  the  extent  of  damage  in  this  case  would  be  5%. 
Therefore,  the  statistics,  such  as  the  mean  and  upper  bound 
limit  for  the  extent  of  damage  provides  a  physically  mean¬ 
ingful  description  of  the  expected  state  of  structural  damage 
due  to  fatigue  cracking  at  any  service  time.  This  quantita¬ 
tive  type  of  information  provides  a  sound  basis  for  evaluat¬ 
ing  structural  durability  requirements  and  for  assessing  dur¬ 
ability  design  tradeoffs  for  metallic  durability-critical 
components . 


SECTION  V 

DEMONSTRATION  OF  DURABILITY  ANALYSIS  METHODS 

The  two-segment  deterministic  stochastic  crack  growth 
approach  (DCGA-SCGA)  for  durability  analysis  described  in 
Section  III  and  documented  in  Volume  I  [2]  is  demonstrated 
in  this  section.  Durability  analysis  methods  for  predicting 
the  crack  exceedance  probability,  p(i,T*),  and  the  cumula¬ 
tive  distribution  of  service  time  to  reach  any  crack  size, 

FT(x  ) (fc) t  are  demonstrated  using:  (1)  coupon  specimen  and 
(2)  the  F-16  lower  wing  skins. 

5.1  DEMONSTRATION  FOR  DOG-BONE  SPECIMENS 

The  advanced  durability  analysis  method  described  in 
Section  III  is  demonstrated  for  both  countersunk  and  straight- 
bore  fastener  holes  in  the  following.  The  initial  fatigue 
quality  is  established  based  on  fractographic  results  obtain¬ 
ed  using  narrow  specimens.  Then,  predictions  for  the  crack 
exceedance  probability,  p(i,*T),  and  cumulative  distribution 
of  service  time  to  reach  a  specific  crack  size,  FT^X  j (t) ,  in 
the  large  crack  size  region  are  made  using  the  DCGA-SCGA. 
Predictions  are  correlated  with  actual  fractographic  results 
obtained  using  wide  dog-bone  specimens. 

5.1.1  Countersunk  Fastener  Holes 

The  initial  fatigue  quality  of  countersunk  fastener 
holes  will  be  determined  using  the  narrow  width  specimen 
(Fig.  5-1)  test  results,  i.e.,  AFXLR4,  AFXMR4  and  AFXHR4  data 
sets.  Then,  the  durability  analysis  prediction  will  be  made 
for  the  test  results  of  wide  width  specimens  (Fig.  5-2) , 
i.e.,  WAFXMR4  and  WAFXHR4  data  sets  where  large  fatigue 
cracks  exist.  Correlations  between  the  theoretical 
predictions  and  test  results  will  be  made  to  demonstrate  the 
validity  of  the  durability  analysis  methodology.  The  narrow 
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Figure  5-1.  N>rro %r  ISt  Bolt  Load  Transfer  Specimen  Deelgn 
(W  ■  i.S**). 


(W  -  1.50")  and  wide  (W  **  3.00")  width  data  sets  used  are 
described  in  Tables  5-1  and  5-2,  respectively. 

The  procedures  used  in  the  demonstration  are  given  as 
follows: 

1.  Use  the  Weibull  compatible  distribution  function  and 
the  pooled  data  sets  (i.e.,  AFXLR4 ,  AFXMR4  and  AFXHR4)  in 
Section  3.3.1  of  Volume  II  [3]  to  define  the  IFQ  of  counter¬ 
sunk  fastener  holes  in  7475-T7351  aluminum.  The  following 
EIFSD  parameters  are  obtained:  xu  ■  0.03", cC*  1.716  and  p  - 
6.308  (see  Table  5-3) . 

2.  Determine  the  crack  growth  rate  parameter  «1  for 
HAFXHR4  and  WAFXHR4  data  sets  in  the  small  crack  size  region, 
(i.e.,  AL-AU  -  0.01"  -  0.05"),  using  the  pooled  Q  values  from 
AFXMR4  and  AFXHR4  data  sets,  respectively.  Determine  the 
crack  growth  rate  parameter  Q2  and  the  corresponding  standard 
deviation  c £  for  WAFXMR4  and  WAFXHR4  data  sets  in  the  large 
crack  size  region  (i.e.,  aQ-AU*  “  0.05"-0.5")  using  the 
fractographic  results  of  WAFXMR4  and  WAFXHR4  data  sets,  res¬ 
pectively. 


3.  Use  the  DCGA-SCGA  to  predict  the  c rack  exceedance 
probability  p(t,T)  in  the  large  crack  size  region  and  the 
distribution  of  service  time  F^^  j(t)  to  reach  a  specific 
large  crack  size  x^. 


4.  Correlate  analytical  predictions  with  the  actual 
test  results  for  two  wide  specimen  data  sets;  i.e.,  WAFVMR4 
and  WAFXHR4. 


WAFXMR4  and  WAFXHR4  data  sets  were  tested  using  the  F-16 
400  hour  spectrum  with  a  maximum  peak  gross  stress  of  34  ksi 
and  40.8  ksi,  respectively.  The  "F-16  400  hour  spectrum"  has 
been  used  extensively  in  recent  years  in  General  Dynamics' 
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Table  5-1.  Description  of  Fnctasrmphic  Data  Sets  Deed  to 

Determine  the  ZFQ  for  Countersunk  Fastener  Holes, 


DM.1  SET 
(Raf.  25) 

SPEC  I)OS 
USED  (4) 

<r(J> 

(tad) 

a 

(%) 

H1DTH 

(la.) 

t 

(la.) 

FASTENER  (2) 

LOAD 

SPECTRUM 

AfSJM 

10/11  <$} 

32 

ii 

L, 

5 

.11 

72 

BM35 

3-#§ 

MS  4 

1 

oo  at 

mm 

vs 

34 

mm 

10/1# 

3# 

\ 

1 

Not*:  (1)  arterial:  7475-77351  altalm 

(2)  Iliad  pell-ttrcegti  rlwt  (cortxnk  hast) 

(3)  Gross  section  stress 

(4)  a/jj  *  So.  of  eertaea  osady total  so.  at  apsiua*  is  feta  art 

(5)  Otirtsd  exact  no.  I  fret  data  at 


Table  5-2 .  Description  of  M&VXMR4  and  *A ROOM  FractoffrapiUc 
Data  Sets. 


DATA  8KT 

Basal 

190. 

CRACKS 

mm 

»J7TTTa(| 

WIDTH 

(In.) 

NlfXNR4 

15 

14 

34 

WAJXW4 

15 

13 

40. S 

Notes t  1.  7475-T73SI  Aluminum 

2.  Ref.  Fig  5-2  for  specimen  design  details. 

3.  Ref.  39 
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IRAD  and  CRAD  research  programs.  However,  this  spectrum 
doesn't  apply  to  P-16  production  aircraft.  Theoretically, 
there  is  no  significant  difference  in  the  peak  stress  at  the 
edge  of  the  fastener  hole  for  narrow  (W  »  1.5")  or  vide  (w  * 
3.0")  specimen  subjected  to  the  same  gross  section  stress. 
The  narrow  specimen  has  a  slightly  larger  net  section  stress 
than  the  wide  specimen.  However,  the  narrow  specimen  has  a 
smaller  stress  concentration  factor  than  the  wide  specimen. 
These  compensating  factors  are  the  reason  the  maximum  peak 
stress  at  the  edge  of  the  fastener  hole  is  virtually  the  same 
for  both  narrow  and  wide  specimens  subjected  to  the  same 
gross  section  stress. 

5. 1.1.1  Estimation  of  Service  Crack  Growth  Parameters. 

EXFSD  parameters  for  countersunk  fastener  holes  based 
on  specimen  data  sets  AFXLR4,  AFXMR4 ,  and  AFXHR4,  are  shown 
in  Table  5-3.  Pooled  Q  values  for  each  of  these  data  sets 
are  also  shown  in  Table  5-3. 

The  crack  growth  rate  parameters  Q1  and  Q2  vary  with  re¬ 
spect  to  service  loading  conditions.  However,  when  all  ser¬ 
vice  loading  conditions  are  identical,  such  as  loading  spec¬ 
tra,  percentage  of  load  transfer,  type  of  fastener  holes, 
etc. ,  except  the  maximum  gross  section  stress  level  <3T,  a  very 
reasonable  model  relating  the  crack  growth  rate  parameter  Q 
and  the  maximum  gross  section  stress  is  given  in  Eg.  (3-32). 

Thus,  if  fractographic  data  sets  are  available  under 
several  different  gross  stress  levels, Cg,  the  empirical  par¬ 
ameters  C  and  V  in  Eq. (3-32)  can  be  determined  using  the  least 
square  fit  procedure.  Then,  the  crack  growth  rate  parameters 
Q1  and  Q2  under  different  gross  stress  levels  can  be  computed 
from  Eq.  (3-32).  For  demonstrative  purpose,  since  applicable 
fractographic  results  in  the  small  crack  size  region  are  av¬ 
ailable  for  AFXLR4,  AFXMR4  and  AFXHR4  narrow  specimen  data 
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sets,  Eq,  3-32  is  used  to  determine  the  crack  growth  rate  pa¬ 
rameters  Q.  for  WAFXMR4  and  WAFXHR4  data  sets  as  well  as  var¬ 
ious  stress  regions  in  the  lower  wing  skin  of  the  F-16  aircraft. 

In  the  small  crack  size  region  of  AL-AU  *  0.01'*  -  0.05”, 

Q  values  versus  gross  stresses  for  the  AFXLR4,  AFXMR4  and 
AFXHR4  data  sets  shown  in  Table  5-3  are  plotted  in  Fig.  5-3 
as  solid  circles.  Using  the  model  of  Eg.  (3-32)  and  a  least- 
squares-fit  procedure,  a  straight  line  is  obtained  in  Fig. 
5-3?  with  C  »  4.829X10"4  and  V  -  6.38.  With  the  values  of  C 
and  V  given  above  as  well  as  the  gross  stresses  for  WAFXMR4 
and  WAFXHR4  data  set,  Q.  values  for  these  two  data  sets  are 

_4  -4 

computed  from  Eq.  £3-32 )as  2.851x10  per  hour  and  9.126x10 
per  hour,  respectively. 

Fractographic  results  available  in  the  large  crack  size 
range,  i.e.,  a0«AU  -  0.05"  -  0.5",  for  AFXLR4,  AFXMR4  and 
AFXHR4  data  sets  are  not  sufficient  to  determine  the  respec¬ 
tive  pooled  Q2  values,  because  the  specimens  for  these  data 
seta  are  only  1.5"  wide.  As  a  result,  the  crack  growth  rate 
parameters  Q2  and  the  corresponding  standard  deviation  <rz 
for  egment  2,  i.e.,  a0-AU* «■  0.05"  -  0.5",  for  WAFXMR4  and 

WAFXHR4  were  determined  using  the  fractographic  results  of 
these  two  data  sets.  Q2  and  °z  values  for  WAFXMR4  and 
WAFXHR4  are  summarized  in  Table  5-4  in  which  the  value  of  Q2 
is  denoted  as  Q. 

5. 1.1. 2  Theoretical/Experimental  Correlations .  Theore¬ 
tical  predictions  for  the  probability  of  crack  exceedance 
p(i,  T") ,  and  the  cumulative  distribution  of  time  to  reach  a 
given  crack  size  FT(X  ^  (t)  ,  for  the  WAFXMR4  and  WAFXHR4  data 
sets  have  been  computed  using  the  DCGA-SCGA.  All  results  are 
based  on  the  following  EIFSD  parameters  for  the  Weibull  com¬ 
patible  distribution  function:  xu*0.03",  o£  «  1.716,  p  * 

6.308  (see  Table  5-3). 
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CRACK  GROWTH  RATE  PARAMETER:  Q  *10 


10 


Figure  5-3.  Crack  Growth  Rate  Parameter  Q  Versus  Gross  Stress 

for  Narrow  Specimen  Data  Sets  (AFXLR4.  AFXMR4 .  AFXHR4). 
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Table  5-4.  Summery  of  Q  and  for  WAFXMR4  and  WAFXHR4 

Data  Seta 


DATA 

% 

n 

WIDTH 

a  -  AU' 

0  KlO4 

(2) 

SET  (1) 

LT 

STRESS 

(ksi) 

(In.) 

0 

(1/Hr.) 

ESI 

15 

14 

34 

M 

SB 

2.906 

3.449 

HA75C4 

15 

13 

40.8 

KH 

3.854 

3.322 

Notts:  (1)  Bef.  Fig.  5-2  for  speciaaa  design  details  (7475-T7351 
aluiima) 

(2}  Bef.  Eg. (3-36) (Natural  log  bvie) 
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The  predicted  probability  of  crack  exceedance  at  7“  m 
11,608  flight  hours  for  WAFXMR4  is  displayed  in  Fig.  5-4  as  a 
solid  curve.  Also  shown  in  this  figure  as  solid  circles  are 
the  actual  test  data  for  comparison.  Further,  the  predicted 
probability  of  crack  exceedance  at  T  *  7,000  flight  hours  for 
WAFXHR4  is  shown  in  Fig.  5-5  as  a  solid  curve  and  the  solid 
circles  denote  the  actual  fractographic  test  results. 

The  predicted  cumulative  distribution  of  service  time  to 
reach  a  crack  size  of  0.73"  for  the  WAFXMR4  data  set  is  dis¬ 
played  in  Fig.  5-6  as  a  solid  curve.  The  actual  fractograph¬ 
ic  results  are  shown  in  the  same  figure  as  solid  circles  for 
comparison.  Similarly,  the  prediction  for  the  cumulative 
distribution  of  service  time  to  reach  a  crack  size  of  0.59" 
for  WAFXHR4  is  shown  in  Fig.  5-7  as  a  solid  curve.  The  solid 
circles  depicted  in  the  same  figure  are  the  actual  fracto¬ 
graphic  test  data  for  comparison.  It  is  observed  from  Figs. 
5-4  to  5-7  that  the  correlations  for  countersunk  fastener 
holes  between  the  experimental  results  and  the  durability 
analysis  predictions  are  very  reasonable. 

5.1.2  Straight-Bore  Fastener  Holes 

The  DCGA-SCGA  for  durability  analysis  is  demonstrated 
for  straight-bore  clearance-fit  fastener  holes  in  7475-T7357 
aluminum  in  this  section,  procedures  for  the  demonstration 
are  given  as  follows. 

1.  The  IFQ  for  straight-bore  clearance-fit  fastener 
holes  is  based  on  the  Weibul  1-compatible  EIF.eD.  Two  narrow 
width  (W  ■  1.5")  specimen  data  sets  (WPF  and  XWPF;  see  Figs. 
5-8  and  5-9,  and  Table  5-5),  a  data  pooling  procedure  and  a 
statistical  scaling  technique  [2]  were  used  to  estimate  the 
EIFSD  parameters  «  and  ^  for  xu  ■  .03".  Results  are  summar¬ 
ized  in  Table  5-6. 
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WBBU1  COMPATBUE 
TAU  1 1.608  HOURS 
ALPHA  1.716PHAI  S3C0XU  Q03 
STD.  Q449.Q  .000288 1.Q1  0.0002906 
2-S6GM04T  DCGA-SCGA  APPROACH 
SCAUNG  FACTOR  4.  POOL©  DATA 


WAPXMR4 


? 

8 

£ 


0.2  h 


Figure  5-4. 


0.35  0.70 

CRACK  SIZE  INCH 


1.05 


9*****a  Pr^dlctad  Crack  R^raadaneit 
at  T  •  11608  Flight  Hour*  for 
D*ta  |a£  Sag  Actual  Photographic  Raauita. 


0.2  r 


WAFXHR4 
WEBULL  COMPATBLE 
TAU  7000  HOURS 
ALPHA  1.710.PHAI  6308.XU  0.03 
STD.  0.322.Q  0009126.Q1  00003854 
2 -SEGMENT  DCGA-SCGA  APPROACH 
SCALING  FACTOR  4.  POOLED  DATA 


0.25  0.50 

CRACK  SIZE  INCH 


0.75 


Figure  5-5.  Corralation  Batwaan  Fradictad  Crack  Escaadanca 
Probability  >t  T  «  7000  Flight  Hours  for 

WAFAHR4  Data  gat  and  Actual  Fractoqraphic  Raaults. 
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CUMULATIVE  DISTIRBUTION 


WAFXA*I4 


8  °-6 


fa  0.4 
0 


*  WEBULL  COMPATBLE 
'  CRACK  SIZE  0.73  N 
ALPHA  1.716PHAI  &3Q8.XU  003 
STD.  0.449,0  .0002861,01  0.0002906 
2-SEGMENT  DCGA-SCQA  APPROACH 
SCALING  FACTOR  4.  POOLED  DATA 


TME,  1000  HOURS 


Figure  5-6. 


RTJT  H 


WAFXHR4 


e 

!  !  “ 

D  H 

“  m  o.« 

M  0 


•  WEBULL  COMPATBLE 
CRACK  SIZE  0.59  N 
ALPHA  1.716.PHAI  6308.XU  0.03 
STD.  0.32 2, Q  0009126.Q1  0.0003854 
2 -SEGMENT  DCGA-SCQA  APPROACH 
SCALING  FACTOR  4.  POOLED  DATA 


Figure  5-7 


TME.  1000  HOURS 


<?9rrtlaUQB  Batwf  a-Prt4lg£t4 

Service  Tlae  to  Reach  0.59”  CrecK  Sine  for  WAF3MR4 
Pete  Set  end  Actual  Frectogrephlc  Results. 
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(«.*..  XWPf.  LYVJPF,  HYWPr ) 


Table  5-5.  Description  of  Fractographic  Data  flats  Used  to 

Determine  tha  IFQ  for  Straight-flora  Fastener  Holes. 


No.  of 

Soeoiaens  Used 

nil 

m 

m 

Fastener 

Load  L 

Spectrua  1 

WPF 

31/33  (2) 

34 

B 

SB 

123^^3 

on 

XWPF 

31/33  (3) 

34 

15 

1.5 

Notes: 


(1)  7475-T7331  Alualnua 

(2)  Delated  fatigue  oraoka  i 2  and  6 

(3)  Deleted  fatigue  oraoka  111  and  16 

(4)  Gross  section  stress  for  peak  spectrua  load 

(5)  Ref.  FHQ  program  [37] 


Table  5-6.  Summary  of  IFQ  Parameters  for  Pooled  Straight-Bore 

Hole  Data  Sets  Based  on  Weibull  Compatible  Distribution 
Function 


DATA  SET 
(1) 

NO. 

SPECIMENS 

AL  -  AU 

Q  xlO4 
(1/Hr. ) 

*u 

0 1 

4> 

l 

jWPF  \ 
|XWPF J 

f 31/337 
\31/33J 

a01"-a05" 

(2.3297 
|3 . 671 J 

q  03M 

4.782 

4.658 

1 

B 

(1) 

(2) 


Notes: 


Ref.  37 

CLSSA  and  "EIFS  fit”  used 


2.  The  crack  growth  rate  model  of  Eqs. (3-16) and (3-17) 
(with  b1  «  b2  ■  1)  and  fractographic  data  for  the  WWPF  data 
set  are  uaed  to  estimate  the  crack  growth  parameter  Q1  and  Q2 
respectively.  Specimen  design  details  for  the  WWPF  data  set, 
shown  in  Fig.  5-10,  are  the  same  as  the  WPF  data  set  in  the 
test  section,  except  that  the  WWPF  specimen  is  wider  (i.e., 
3.0"  width).  Such  specimens  are  wide  enough  to  provide  frac¬ 
tographic  data  in  the  large  crack  size  region.  Specimens  for 
the  WWPF  data  set  were  fatigue  tested  to  failure  using  the 
same  load  spectrum  (F-16  400  hour)  and  maximum  peak  (gross) 
stress  level  (i.e.,  34  ksi)  as  the  "WPF"  data  set.  In  the 
present  demonstrations,  AL-AU  ■  0.01"  -  0.05 "is  used  for  the 
small  crack  size  region  (i.e.,  "segment  1")  and  ag-AU* -0.05" 
-  1"  is  used  for  the  large  crack  size  region  (i.e.,  "segment 
2").  Results  for  Q^,  Q2  and  c Tz  f or  the  WWPF  data  set  are 
summarized  in  Table  5-7. 

3.  Theoretical  predictions  for  the  probability  of  crack 
exceedance,  p(i,T'),  at  service  time  T"  18,400flight  hours, 
are  shown  in  Fig.  5-11  for  the  DCGA-SCGA.  Experimental  re¬ 
sults  are  denoted  as  solid  circles  for  comparison. 

4.  Theoretical  predictions  for  the  cumulative  distri¬ 
bution  of  service  time  to  reach  a  crack  size  x^^  ■  are 
shown  in  Fig.  5-12  for  the  DCGA-SCGA.  Experimental  results 
for  the  WWPF  data  set  are  plotted  as  plus  signs  (+)  for  com¬ 
parison. 

The  theoretical  predictions  shown  in  Figs.  5-11  and  5-12 
correlate  well  with  actual  test  results  for  both  the  small 
and  large  crack  size  regions.  Hence,  the  DCGA-SCGA  for  dur¬ 
ability  analysis  can  be  used  to  assess  functional  impairment 
due  to  excessive  cracking,  fuel  leaks  and/or  ligament  break¬ 
age. 
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Table  5-7 .  Summary  of  Q  and  Values  for  WWPP  Data  Sat 


t- 


I 


A 


WWPF 

WBBULL  COMPATBLE 
TAU  18.400  HOURS.  AO  Q06  NCH 
ALPHA  4782.  PHA1  46881  XU  003 
STD  0177.  Q1  .000274ZQ2  .0003124 
2-SK»»IT  OCGA-SCGA  APPROACH 
SCAUNG  FACTOR  1.  XWPP  4 


078 


CRACK  SCE 


Figure  5-11.  Correlation  Between  Predicted  Crack  gxceedanc* 
Probability  Baaei  on  the  DCGA-SCGA  atT  ■  18,400 
Flight  Hour*  for  WWPF  Data  Set  and  Actual  Fract- 
ograpblc  Reaulta. 


* 
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Figure  5-12.  Correlation*  Between  Theoretical  Predictions  and 
Experimental  Results  (WWPF  Data  Set)  for  Cumul¬ 
ative  Diatributlon  of  Service  Time  to  Reach  Cracft 
5Ise  a,  « o.  Su  Based  onft&SA-^CGA. 


v 
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5.2  DEMONSTRATION  FOR  THE  F-16  LOWER  WING  SKINS 


A  durability  analysis  of  tha  F-16  lower  win?  skins  has 
been  previously  reported  [1,17,20-22].  This  analysis  was 
concerned  with  relatively  small  fatigue  cracks  (e.g.,  < 

for  excessive  cracking  and  reflected  the  one-segment 
DCGA  [1,16]. 

A  durability  analysis  of  the  F-16  lower  wing  skin  for 
functional  impairment  is  conducted  herein  using  the  two-seg¬ 
ment  DCGA-SCGA.  The  two-segment  DCGA-SCGA  is  demonstrated 
and  evaluated  in  the  following  and  in  Volume  II  [3].  Predic¬ 
tions  will  be  correlated  with  results  from  the  F-16  wing  dur¬ 
ability  test  article.  The  F-16  wing  box  assembly  is  shown  in 
Fig.  5-13  and  stress  regions  for  the  lower  wing  skin  are 
shown  in  Fig.  5-14. 

A  full-scale  F-16  wing  durability  test  was  conducted  us¬ 
ing  the  F-16  1000  hour  spectrum,  consisting  of  two  500  hour 
blocks.  After  fatigue  testing  to  16,000  flight  hours,  a 
tear-down  inspection  was  performed.  All  fastener  holes  in 
both  lower  wing  skins  (i.e.,  3228  holes)  were  inspected  using 
the  eddy  current  technique.  Each  fastener  hole  with  a  crack 
indication  was  broken  open  and  a  fractographic  analysis  was 
performed.  Tear-down  inspection  and  fractographic  results 
are  documented  in  Ref.  38. 

The  following  procedures  are  used  to  demonstrate  and  ev¬ 
aluate  the  two-segment  DCGA-SCGA  using  the  F-16  lower  wing 
skins  for  /  the  durability  analysis  for  functional  impairment 
associated  with  large  through-the-thickness  cracks. 

1.  The  IFQ  is  based  on  the  fractographic  results  from 
AFXLR4,  AFXMR4  and  AFXHR4  data  sets.  The  EIFSD  parameters, 
based  on  the  fractographic  data  in  the  small  crack  size  range 
AL-AU  -  0.01”  -  0.05”  andJt-  4  for  each  of  the  three  data 
sets  heve  been  obtained  in  the  previous  exaiu;*le  to  define  the 
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IFQ  of  countersunk  fastener  holes;  with  the  results  xu  « 
0.03",  1.716  and  f  -  6.308,  see  Table  5-3. 

2.  The  F-16  lover  wing  skin  is  divided  into  10  stress 
regions  as  shown  in  Fig.  5-14.  The  stress  level  and  the  num¬ 
ber  of  fastener  holes  in  each  stress  region  are  shown  in 
Table  5-8. 


3.  The  crack  growth  rate  parameter,  Q,  for  segment  1, 
in  each  stress  region  are  determined  using:  (i)  the  avail¬ 
able  pooled  Q  values  from  the  AFXLR4,  AFXMR4  and  AFXHR4  data 
sets  (see  Table  5-2;  AL-AU  -  0 . 01"-0 . 05" ) ,  and  (ii)  the  model 
for  Q  as  a  function  of  stress  given  by  Eq. (3-32 Results  of 
the  model  parameters  C  and  V  in  Eq. (3-32)  obtained  from  three 
data  sets  (AFXLR4 ,  AFXMR4  and  AFXHR4)  have  been  computed  in 
the  previous  example,  Fig.  5-3,  and  they  are  shown  in  Fig.  5- 
15,  Frame  A. 

4.  The  crack  growth  rate  parameters,  Q2,  for  Segment  2 
in  each  stress  region  are  determined  using  available  wide 
specimen  fraotographic  results  from  WAFXMR4  and  WAFXHR4  data 
sets  in  a0-AU#  -  .05"  -  .5"  along  with  Eq.  (3-32).  The  model 
parameters  c  and  V  obtained  from  WAFXMR4  and  WAFXHR4  data 
sets  are  shown  in  Fig  5-15,  Frame  B. 

5.  Predictions  for  p(i ,f)  in  each  stress  region,  based 
on  the  two-segment  DCGA-SCGA,  are  computed  using  Eqs.(3-19), 
(3-23 )  ,(3-29)  and  (3-30)  . 

6.  From  the  predicted  crack  exceedance  probability, 
p(i, T)  and  the  number  of  fastener  holes  in  each  stress  re¬ 
gion,  the  statistics  for  the  number  of  cracks  exceeding  some 
crack  sizes  in  the  entire  lover  wing  skin  are  computed  using 
the  Binomial  distribution  Eqs.  (3-37) to  (3— 40 )  [e.g. ,  40]. 

7.  Theoretical  predictions  are  correlated  with  actual 
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Table  5-8.  Stress  Levels  end  Nuaber  o£  Fastener 
Holes  for  F-16  Lower  Wing  Skin 


MAX.  STRESS 
LEVEL  (ksi) 

NO.  OF  FASTENER 
HOLES 

i 

28.3 

59 

IX 

27.0 

320 

XXX 

24.3 

680 

IV 

16.7 

469 

V 

28.4 

8 

VX 

29.2 

30 

VII 

32.4 

8 

VIIX 

26.2 

6 

XX 

26.2 

12 

X 

25.7 

20 

1614 
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an-AU  -  .05*'-. 5" 


*  jrfuAnuim 
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/  |  °i  • c<rV 

« 

-6 

'  C  -  1.234*10 

|  V  -  1.549 

1  I 

For  each  of  10  stress  regions 


Figure  5-15.  General  Approach  for  Estimating  Service  Creek 
Growth  Parueeters  Q  ^  andi  Q2  . 
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test  results  from  the  F-16  durability  test  article.  Results 
will  be  plotted  in  a  useful  format  for  evaluating  the  two- 
segment  DCGA-DCGA  and  the  DCGA-SCGA  for  durability  analysis. 

The  same  three  fractographic  data  sets,  i.e.  AFXLR4, 
AFXMR4  and  AFXHR4,  were  used  to  determine  the  EIFSD  paramet¬ 
ers  in  the  previous  [1,17,20-22]  and  present  durability  an¬ 
alyses  for  F-16  lower  wing  skin.  However,  different  oe  and 
$  values  for  xu  -  0.03"  are  obtained  in  the  present  analysis 
due  to  the  difference  in  the  following:  (1)  fractographic 

crack  size  AL-AU  ranges  used  for  determining  EIFSD  and  (2) 
fractographic  data  processing  methods/screening  considera¬ 
tions  used.  The  resulting  EIFSD  parameter  values  are  xu  - 
0.03",  oL  -  1.716  and  6.308  (see  Table  5-3). 

In  the  previous  durability  analysis  [1],  terminal  crack 
size  dimensions  in  fastener  holes  were  based  on  initial  meas¬ 
urements  of  the  fracture.  In  the  present  durability  analy¬ 
sis,  however,  terminal  crack  sizes  ware  based  on  the  fracto- 
graphy.  The  final  crack  dimension  based  on  the  fractography 
are  more  accurate  than  the  initial  fracture  surface  measure¬ 
ments.  There  are  small  differences  between  the  initial  crack 
size  dimensions  and  those  based  on  the  fractography.  As  a 
result  of  these  differences,  the  experimental  results  for  the 
average  number  of  fastener  holes/skin  (for  both  wing  skins) 
with  a  crack  size  >  is  14.5  holes  (fractography)  versus 

16.5  holes  (initial  measurements). 

The  F-16  durability  test  article  was  fatigue  tested  to 
flight  hours  using  the  F-16  1000  hour  load  spectrum. 
This  preliminary  spectrum  included  two  500-hour  blocks.  The 
F-16  400  hour  loading  spectrum  has  been  used  extensively  in 

recent  years  for  General  Dynamics  IRAD  and  CRAD  research  pro¬ 
grams.  This  spectrum  is  slightly  more  severe  than  the  F-16 
hour  spectrum  but  it  doesn’t  apply  to  F-16  production 
aircraft.  It  is  assumed  for  durability  analysis  purposes 
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that  the  coupon  fractographic  results  (i.e.,  AFXLR4 ,  AFXMR4, 
AFXHR4 ,  WAFXMR4  and  WAFXHR4)  based  on  the  F-16  400  hour  spec¬ 
trum  can  be  applied  for  the  prediction  of  the  F-16  durability 
test  article. 

The  F-16  lower  wing  skins  contain  several  cutouts.  How¬ 
ever,  the  present  durability  analysis/correlation  covers  only 
fatigue  cracks  in  fastener  holes. 

5.2.1  Estimation  of  Service  Crack  Growth  Parameters 

The  service  crack  growth  parameters  Q^,  Q2  and  d*z  were 
estimated  for  the  small  (i.e.,  AL-AU  ■  0.01"  -  0.05") and  large 
crack  size  region  (i.e.,  a0-Au‘  «  0.05"-0.5")  for  each  of 
the  ten  stress  regions.  A  general  approach  for  estimating  Q1 
and  Q2  is  described  in  Fig.  5-15.  In  the  small  crack  size 
region,  Q1  values  for  the  AFXLR4 ,  AFXMR4  and  AFXHR4  data  sets 
were  obtained  previously,  Table  5-3  and  Fig.  5-3.  From  these 
Q1  values,  the  constants  C  and  V  in  Eq. ( 3-32) were  determined 
using  a  least-squares  fit  procedure  (Fig.  5-3) .  Then,  Q1 
values  in  each  of  the  ten  stress  regions  are  computed  from 
Eq.  3-32,  and  the  results  are  shown  in  Table  5-9. 

A  similar  approach  to  that  described  above  was  used  for 
the  large  crack  size  region  to  estimate  Q2  for  each  of  the 
ten  stress  regions.  In  this  case,  fractographic  results  of 
the  WAFXMR4  and  WAFXHR4  data  sets  (see  Table  5-4)  were  used 
to  estimate  the  constants  C  and  V  in  Eq. (3-32).  Result#  are 
shown  in  Table  5-9  and  in  Fig.  5-16. 

In  practice,  suitable  fractographic  data  may  not  be  av¬ 
ailable  to  estimate  Qx  and  Q2.  In  such  cases,  an  analytical 
crack  growth  program  [e.g.,  31,32]  can  be  used  to  estimate 
the  crack  size  versus  time  information  needed  to  establish  Q1 
and  Q2  for  given  durability  analysis  conditions  (e.g.,  stress 
level,  load  spectrum,  %  bolt  load  transfer,  etc.).  Refer  to 


5-27 


Table  5-9.  Summary  of  Crack  Growth  Rata  Parameters 
for  Each  Stress  Region. 


STRESS 

REGION 


MAX.  STRESS 

NO.  07  FASTENER 

Q,  *104(l) 

LEVEL  (ksi ) 

HOLES 

(1/HR. ) 

28.3 

59 

.884 

27.0 

320 

.655 

24.3 

680 

.334 

16.7 

469 

.030 

26.4 

8 

.904 

29.2 

30 

1.080 

32.4 

8 

2.097 

26.2 

8 

.541 

26.2 

12 

.541 

25.7 

20 

.478 

2.187 

2.033 

1.727 

.966 

2.199 

2.296 

2.697 

1.941 

1.941 

1.884 


Notes: 


(1)  Segment  1:  AL-AU  >0.01"  -  0.5" 
cx  >  4.829*10’14;  V,  >  6.380 

(2)  Segment  2:  AL-AU  *  7).  05"  -  0.5" 
C2  ■1.234*10'"®  i  V2  -  1.549 


trail  STRCSf  .  Mt»a 


a  0  -  Ml  =  .05"  -  .5" 


«AfXW4 


WnMM 


Qi  ■  1.234*io'tr 1,349 
C  Sails  (ksl) 


to  10  *0  90 

attest  itrcss  um 

Figure  5-16.  Crack  Growth  Rate  Parameter  Q  Versus 
Gross  Stress  for  Wide  Specimen  Data 
Sets  (WKFXMR4  and  WAF3blR4). 


Section  III  herein  and  to  Vols.  I  [2]  and  II  [3]  for  further 
details. 


5.2.2  Theoretical/Experimental  Correlations 

Probability  of  crack  exceedance  predictions  p(i, T)  at 
X «  16, 000  flight  hours  for  five  different  crack  sizes  (i.e., 
x1  »  0.03”,  0.05",  0.1",  0.2”  and  0.3")  are  shown  in  Table  5- 
10  for  the  two-segment  DCGA-SCGA.  The  average  number  of 
fastener  holes  in  each  stress  region,  T?(i,T)  with  a  crack 
size  greater  than  x^  at  f  -  16,000  flight  hours  is  also 
shown  in  this  Table.  The  analysis  for  the  DCGA-SCGA  was  con¬ 
ducted  using  (C  ■  0,3  (natural  log  basis),  which  is  quite 
reasonable  for  countersunk  fastener  holes  in  7475-T7351  alum¬ 
inum  [3,6]. 

Predictions  for  the  average  number  of  fastener  holes  in 
the  entire  lower  wing  skin  with  a  crack  size  >  x.  at  16, 000 
flight  hours,  E( T-) ,  and  the  standard  deviation,  d^CT) ,  are 
shown  in  Table  5-11  for  both  the  DCGA-DCGA  {see  Vol.  II  [3]) 
and  the  DCGA-SCGA.  Z{T)  and  0[(T)  values  are  computed 
based  on  the  Binomial  distribution,  Eqs. (3-39)  and  0-40).  The 
tear-down  inspection  results  based  on  the  average  of  two  low¬ 
er  wing  skins  are  shown  in  the  same  table  for  comparison. 

Theoretical  predictions  for  the  average  number  of  fas¬ 
tener  holes,  L(T) ,  with  a  crack  size  >  x1  atT"-  flight 

hours  in  the  entire  lower  wing  skin  are  plotted  in  Fig.  5-17 
for  both  of  the  two-segment  crack  growth  approaches  for  com¬ 
parison  purposes.  In  this  figure,  the  results  for  the  DCGA- 
DCGA  and  the  DCGA-SCGA  are  depicted  by  a  solid  curve  and  a 
dashed  curve,  respectively.  Results  for  both  approaches  are 
identical  for  the  crack  size  xx  <  0.05"  in  the  first  crack 
growth  segment.  The  tear-down  inspection  results  are  shown 
in  Fig.  5-17  as  solid  circles  for  comparison.  These  results 
reflect  the  average  extent  of  damage  for  a  lower  wing  skin 
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Table  5-10.  Crack  Exceedance  Probability  and  Average  Nuo&er 
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Table  5-11.  statistic*  for  Nun be r  of  Fastener  Holaa 
with  Crack  Siaa  Exceeding  *1  i»  F-16 
Lower  Wing  Skin  for  Both  DCGA-DCGA  and 
DCOA-SCGA. 


*1 

(IN.) 

mm 

BXFBKIMfNTAL 
RESULTS  (AVE.) 

mssm 

mtm 

agog 

Kan 

0.03 

35.  SO 

5.800 

5.800 

14.5 

O.OS 

10. SI 

3.185 

10.81 

3.185 

9.5 

0.1 

5.37 

2.258 

5.38 

2.262 

7.0 

1.99 

1.379 

2.19 

1.450 

1.0 

0.3 

1.00 

.977 

1.24 

1.097 

0.5 
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P  *  Exceedance  Probability 


GKACK  S IZX.  INCH 


Figure  5-17. 


Correlations  between  Theoretical  FredJ 
and  Experimental  Result a  for  Flehter 


based  on  the  total  extent  of  damage  for  both  lower  wing  skins 
combined.  Since  the  number  of  details  in  each  stress  region 
is  large,  it  is  reasonable  to  approximate  the  binomial  dis¬ 
tribution  by  the  normal  distribution.  Hence,  the  predicted 
average  extent  of  damage,  L(T)  corresponds  to  an  exceedance 
probability  of  P  -  0.5,  see  Fig.  5-17. 

The  extent  of  damage  estimate  for  an  exceedance  prob¬ 
ability  of  P  -  0.05  is  also  plotted  in  Fig.  5-17  as  a  solid- 

-dashed-solid  curve  ( - ) .  This  curve  represents  the 

estimated  upper  bound  limit  for  the  extent  of  damage  with  an 
exceedance  probability  P  ■  0.05.  It  is  computed  from  L (7")  + 
1.65<£(D  where  L(T')  and d£(T)  values  are  shown  in  Table  5-11 
for  the  DCGA-SCGA. 

To  illustrate  the  usefulness  of  the  extent  of  damage 
concept  consider,  for  example,  the  extent  of  damage  at  x1  - 
0.3"  in  Fig.  5-17.  The  (predicted)  probability  is  50%  (i.e., 
P-0.5)  that  1.24  fastener  holes  will  have  a  crack  size  ex¬ 
ceeding  x1  -  0.3"  ;  whereas,  the  probability  is  5%  (i.e., 
P-0.05)  that  3.05  fastener  holes  will  have  a  crack  size  larg¬ 
er  than  Xj^  «  0.3"  atT-16000  flight  hours.  Therefore,  the 
durability  analysis  provides  quantitative  estimates  of  the 
extent  of  damage  mean  and  upper  bound  limit.  This  type  of 
information  provides  a  physical  description  of  the  state  of 
damage  for  a  durability-critical  component  and  a  logical 
basis  for  estimating  structural  maintenance/ repair  require¬ 
ments  and  costs. 


5.2.3  Discussion  of  Results 

The  two-segment  DCGA-SCGA  has  been  demonstrated  and  ev¬ 
aluated  using  fractographic  results  for  both  coupon  specimens 
and  lower  wing  skins  from  a  fighter  aircraft.  This  approach 
was  evaluated  for  fatigue  cracking  in  both  straight  bore  and 
countersunk  fastener  holes  with  clearance- fit  fasteners.  Re- 
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suits  for  two  different  two-segment  durability  analysis  ap¬ 
proaches  (i.e. ,  DCGA-DCGA  and  OCGA-SCGA)  were  compared  for 
the  lower  wing  skin  demonstration.  Both  approaches  are  con¬ 
sidered  reasonable  for  evaluating  functional  impairment  due 
to  fuel  leakage/ ligament  breakage  in  metallic  aircraft  struc¬ 
tures.  However,  the  DCGA-SCGA  is  recommended  for  durability 
analysis  because  predictions  are  more  accurate  and  slightly 
more  conservative  than  those  based  on  the  DCGA-DCGA.  Exten¬ 
sive  demonstrations  for  the  DCGA-DCGA  were  given  in  Volume  II 
[3]. 


The  stress  level  in  each  stress  region  is  important  for 
crack  growth  predictions.  Therefore,  the  stress  analysis  for 
durability-critical  components  should  reflect  appropriate 
finite  element  grid  sizes  to  obtain  the  stress  analysis  ac¬ 
curacy  desired  for  each  stress  region. 


5-34 


4 


.4 


SECTION  VI 

DURABILITY  ANALYSIS  SOFTWARE 

Software  is  available  for  implementing  the  advanced  dur¬ 
ability  analysis  method  described  in  Section  II  of  this  Vol¬ 
ume  (IV)  and  in  Volume  I  [2].  A  comprehensive  software 
user's  guide  is  given  in  Volume  V  [24]. 

6.1  SOFTWARE  DESCRIPTION 


The  advanced  durability  analysis  software  includes  six 
programs  in  "GWBASIC" .  The  purpose  of  each  program  is  des¬ 
cribed  in  Table  6-1..  All  programs  can  be  implemented  on  an 
IBM  or  IBM-compatible  personal  computer. 

Software  is  available  for  plotting  the  fractographic 
data  for  any  crack  size  or  time  range  and/or  durability  an¬ 
alysis  results  for  FT(X  )  (t) ,  p(i,T)  or  F&(t)(x).  A  plott¬ 
ing  capability  is  available  for  the  following  durability  an¬ 
alysis  options:  (1)  DCGA,  (2)  DCGA-DCGA,  and  (3)  DCGA-SCGA. 
Plots  can  be  obtained  with  or  without  correlating  data.  Ty¬ 
pical  example  plots  are  shown  in  Fig.  6-1. 


6.2  SYSTEM  REQUIREMENTS 

The  advanced  durability  analysis  software  is  programmed 
in  "GWBASIC".  It  runs  on  the  IBM  PC  and  IBM-compatible  sys¬ 
tems  with  the  following  minimum  configuration: 


Memory: 

operating  System: 
Graphics  Monitor: 
Disk  Drive: 
Printer: 

Graphics  Program: 


64 OK  RAM 

MS-DOS  Version  2.0  or  Later 
Monochrome  or  Color 
1  Double  Sided  Disk  Drive 
IBM  or  IBM-Compatible  Graphics 
Printer 

Need  Special  "GRAPHICS "  Program  for 
Doing  Screen  Prints  of  Graphic 
Display 
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Table  6-1.  Description  of  Durability  Analysis 

Software  for  IBM  or  IBM-Compatible  PC. 


M(03MM 

rZLSNMtS 

1 

PURPOSE 

"fract" 

Save  or  read/priat  out  fractographic 
data  on  S  1/4"  floppy  disk 

•’scitsxr* 

study  the  character  and  quality  of  a 
fractographic  data  set  (tabulate  data 
and  plot  fractography) 

MQSZAT** 

coapute  pooled  Q  and  for  a  given 

fractographlo  data  set 

t 

"VCIfC" 

Kstlaate  CXPfO  paraaeters  for  Welbull 
coapa tibia  distribution  function 

"PLOt" 

Plot  fraetographlc  data  and/or  dur¬ 
ability  analysis  results 

"ARAL" 


Make  durability  analysis  predictions 


( ti  /  ei  > iff 


for  Durabllit 
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CONCLUSIONS  AND  RECOMMENDATIONS 

7 . 1  CONCLUSIONS 

1.  A  comprehensive  probabilistic  durability  analysis 
approach  has  bssn  developed  for  metallic  aircraft  structures. 
It  appliss  to  ths  crack  growth  accumulation  in  any  type  of 
structural  detail  (e.g.,  fastener  holes,  cutouts,  fillets, 
etc.).  The  approach  has  been  verified  for  clearance-fit  fas¬ 
tener  holes  in  7475-T7351  aluminum  at  two  levels:  (1)  coupon 
specimens,  and  (2)  full-scale  aircraft  structure.  Very  reas¬ 
onable  durability  analysis  results  have  been  obtained,  in¬ 
cluding  damages  due  to  both  small  cracks  (e.g.,  <  0.05")  and 
large  through-the-thickness  cracks  (e.g.,  >  0.5"). 

2.  It  has  been  shown  that  the  initial  fatigue  quality 
(IFQ)  of  both  straight-bore  and  countersunk  fastener  holes 
with  clearance-fit  fasteners  can  be  reasonably  estimated  us¬ 
ing  fractographic  results  from  coupon  specimens  and  that  the 
IFQ  can  be  represented  by  an  equivalent  initial  flaw  size 
distribution  (EIFSD) .  Furthermore,  it  has  been  demonstrated 
that  the  IFQ  of  fastener  holes  in  full-scale  structures  can 
be  defined  using  coupon  specimens. 

3.  The  probabilistic  durability  analysis  approach  deve¬ 
loped  can  be  used  to  "quantify"  structural  durability  in 
meaningful  terms,  such  as:  (1)  probability  of  crack  exceed¬ 
ance  at  any  service  time,  (2)  probability  of  functional  im¬ 
pairment  at  any  service  time,  (3)  cumulative  distribution  of 
service  time  to  reach  any  given  crack  size,  (4)  extent  of 
damage,  and  (5)  structural  wearout  rate.  Since  the  probabil¬ 
istic  approach  developed  accounts  for  the  fatigue  crack 
growth  accumulation  in  each  structural  detail  susceptible  to 
fatigue  cracking  in  service,  it  is  referred  to  as  a  "qu^uti- 
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tative  durability  analysis  approach."  Tha  extant  of  damage 
prediction  at  a  given  service  time  is  defined  by  the  statis¬ 
tics,  such  as  the  average  and  standard  deviation,  of  the  num¬ 
ber  of  structural  details  expected  to  exceed  functional  im¬ 
pairment  crack  size  limits.  This  quantitative  prediction 
provides  an  effective  basis  for  evaluating  functional  impair¬ 
ment,  economic  life  and  structural  vearout,  and  trade-offs  as 
a  function  of  the  design  and  service  variables. 

4.  The  probabilistic  durability  analysis  approach  is  a 
powerful  "durability  design  tool."  It  gives  the  user  new 
durability  analysis  capabilities  and  features  not  provided  by 
the  existing  deterministic  crack  growth  approach  based  on  the 
"worst  case"  detail  within  a  group  of  details.  The  probabil¬ 
istic  durability  analysis  method  is  not  intended  to  complete¬ 
ly  replace  the  deterministic  crack  growth  approach  in  the 
durability  design  process.  The  deterministic  crack  growth 
approach  will  continue  to  be  a  valuable  tool  for  durability 
analysis  -  primarily  during  the  preliminary  design  process. 
Since  a  deterministic  crack  growth  analysis  provides  informa¬ 
tion  only  for  the  "worst  case"  detail  within  a  group  of  de¬ 
tails,  it  cannot  provide  the  "extent  of  damage"  type  informa¬ 
tion  for  the  entire  population  of  structural  details. 

5.  Equivalent  initial  flaw  sizes  (EIFSs)  are  determined 
by  back-extrapolating  fractographic  results.  Since  the  frac- 
tographic  data  depends  on  the  testing  conditions  (e.g.,  load 
spectrum,  fastener  holes,  cutout,  etc.),  EIFSs  are  not 
strictly  "generic."  However,  EIFSD  parameters  can  be  esti¬ 
mated  for  different  fractographic  data  sets  using  the  data 
pooling  and  statistical  scaling  procedures.  It  has  been  con¬ 
clusively  shown  that  the  EIFSD  based  on  given  fractographic 
data  sets  can  be  used  to  obtain  very  reasonable  durability 
analysis  predictions  for  the  other  data  sets  and  full-scale 
aircraft  structure  for  clearance-fit  fastener  holes  (both 
straight-bore  and  countersunk)  in  7475-T7351  aluminum.  It 
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should  ba  claar  that  an  EIFSD  doas  not  nacaaaarily  contain 
tha  "rogua  flaw." 

6.  Whan  an  EIFSD  la  grown  forward  to  a  aalactad  aarvica 
tine,  tha  aarvica  crack  growth  should  ba  consistent  with  tha 
"basis"  for  tha  EIFSa.  Therefore,  tha  analytical  crack 
growth  program  uaad  [e.g.,  31,32]  should  ba  "tunad"  or  "curva 
fittad"  to  tha  EIFS  master  curves  reflected  in  tha  EIFSD. 

7.  Probabilistic-based  durability  analysis  methods  [2, 
3,5-7]  are  now  sufficiently  developed  and  demonstrated  for 
immediate  applications  to  metallic  airframes.  An  updated 
durability  design  handbook  and  software  for  an  IBM  or  IBM- 
compatible  PC  are  available  for  implementing  the  advanced 
durability  analysis. 

8.  A  "natural  fatigua  crack"  data  base  for  estimating 
the  initial  fatigue  quality  of  structural  details  can  be  ac¬ 
quired  as  a  part  of  tha  Aircraft  Structural  Integrity  Program 
(ASIP)  test  plan.  For  example,  by  not  preflawing  structural 
details  in  test  specimens,  "natural  fatigua  crack"  data  can 
be  obtained — thereby  satisfying  data  requirements  for  both 
durability  and  damage  tolerance.  Additional  testing  and 
fractographic  evaluations,  beyond  the  normal  ASIP  effort,  may 
be  needed  to  define  IFQ,  depending  on  the  desired  confidence 
level  and  circumstances.  IFQ  data  requirements  can  be  read¬ 
ily  incorporated  into  the  ASIP  test  plan  to  minimize  the  cost 
and  time  for  acquiring  the  requisite  data  base. 

9.  The  stress  level  for  each  stress  region  is  important 
for  crack  growth  predictions.  Therefore,  the  stress  analysis 
for  durability-critical  components  should  reflect  appropriate 
finite  element  grid  sizes  to  obtain  the  desired  stress  analy¬ 
sis  accuracy  for  each  stress  region. 

10.  Probabilistic  durability  analysis  methodologies  de- 
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veloped  can  be  extended  to  establish  the  optimal  inspection/ 
repair/replacement/proof  test  maintenance  for  life  management 
of  metallic  aircraft  structure.  The  extension  can  be  made 
based  on  some  fundamental  research  efforts  appearing  in  the 
literature  [e.g,,  18,  35-36,  50-58]. 

7.2  RECOMMENDATIONS 

1.  The  advanced  durability  analysis  method  developed 
under  this  program  should  be  used  for  future  durability  an¬ 
alyses  for  metallic  airframes.  Structural  durability  can  now 
be  quantitatively  accounted  for  in  the  durability  design  pro¬ 
cess. 


2.  Recommendations  for  durability  analysis  are  as  fol¬ 
lows:  (l)  define  the  equivalent  initial  flaw  size  distribu¬ 
tion  (EIFSD)  using  fractographio  data  in  the  small  crack  size 
region  (e.g.,  o.oi"-o.05") ,  (2)  use  fractographio  data  pool¬ 
ing  procedure  and  statistical  scaling  technique  to  estimate 
the  EIFSD  parameters  in  a  "global  sense"  for  a  "single  hole 
population"  basis,  and  (3)  use  the  two-segment  deterministic- 
stochastic  crack  growth  approach  (DCGA-SCGA)  to  predict  the 
extent  of  damage  in  the  entire  durability  critical  component; 
the  two-segment  deterministic  crack  growth  approach  (DCGA- 
DCGA)  is  also  reasonable  but  it  is  slightly  less  conservative 
than  the  DCGA-SCGA. 

3.  The  recommended  changes  in  Air  Force  philosophy  and 
durability  design  requirements  described  in  Volume  IV  [4] 
should  be  adopted.  This  will  allow  the  full  potential  of  the 
probabilistic  durability  analysis  approach  to  be  utilized  in 
the  design  and  analysis  of  future  metallic  aircraft  struc¬ 
tures. 


4.  The  advanced  durability  analysis  approach  developed 
under  this  program  should  be  investigated  for  other  structur- 
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al  detail*  and  considerations.  For  example,  the  life  en¬ 
hancement  effects  of  fastener  hole  cold  working,  interference 
fit  fasteners,  press  fit  bushings,  etc.,  on  initial  fatigue 
quality  should  be  investigated.  Similarly,  the  initial  fa¬ 
tigue  quality  of  structural  details,  such  as  cutouts,  lugs, 
fillets,  etc.,  should  be  investigated.  Suitable  test  speci¬ 
mens  should  be  developed  and  standardized  for  acquiring  init¬ 
ial  fatigue  quality  data  for  those  structural  details  to  be 
included  in  the  durability  analysis. 

5.  Future  ASIP  test  plans  should  be  designed  to  provide 
data  for  initial  fatigue  quality,  durability  and  damage  tol¬ 
erance.  SeleCw^u  fatigue  tests  should  be  conducted  using 
specimens  witMut  intentional  preflaws  so  that  "natural  fa¬ 
tigue  crack"  data  can  be  obtained.  This  approach  should  be 
used  to  minimize  cost  and  time  for  acquiring  the  requisite 
IFQ  data  base. 

6.  The  meaning  and  limitations  of  EIFSs  and  an  EIFSD 
must  be  emphasized.  In  particular,  all  EXFSs  should  be  grown 
forward  consistent  with  the  basis  for  the  EIFSD.  The  EIFSD 
should  not  be  grown  forward  using  an  analytical  crack  growth 
program  without  tuning  and  considering  the  basis  for  the 
EIFS. 


7.  All  aerospace  contractors  should  use  the  same  method 
to  define  EXFSs  for  different  materials  and  structural  de¬ 
tails  so  that  compatible  EXFSs  can  be  obtained.  The  "Qa(t) 
model"  (Eg.  (3-2  )  with  b«l)  is  reasonable  for  determining 
EXFSs.  This  model  or  some  other  suitable  model  should  be 
used  to  standardize  the  way  EXFSs  are  determined.  Then,  for 
a  given  fractographic  data  set,  fractographic  crack  size 
range  (AL-AU)  and  the  same  analysis  procedure,  all  contrac¬ 
tors  will  obtain  the  same  EXFSs.  By  standardizing  the  way 
EXFSs  are  determined,  EIFSs  from  various  sources  can  be  di¬ 
rectly  compared  -  thereby  providing  a  means  for  cataloging 
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and  utilizing  existing  data  from  various  sourcas  to  estimate 
tha  initial  fatigua  quality  of  structural  details. 

8.  Initial  fatigua  quality  and  tha  initial  flaw  size 
distribution  vary  with  respect  to  material,  type  of  fastener 
hole,  structural  details,  manufacturing  processes,  etc.  For 
example,  tha  statistical  dispersion  of  EIFSD  for  countersunk 
holes  is  significantly  larger  than  that  for  the  EIFSD  for 
straight-bora  holes  for  clearance-fit  fasteners  in  the  same 
material  in  which  the  holes  ware  drilled  using  comparable  me¬ 
thods  . 


9.  The  probabilistic  durability  analysis  approach 
should  be  investigated  for  discriminating  "quality"  at  three 
levels:  (1)  material,  (2)  manufactured  detail,  and  (3)  com¬ 
ponent.  Of  particular  interest  is  tha  following  question: 
"How  does  improvement  in  initial  material  quality  translate 
into  improvement  in  life  of  actual  aircraft  components?" 
This  research  can  be  built  on  the  advancements  made  under 
this  program  and  the  work  conducted  by  ALCOA  [e.g.,  48,49]. 
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DEFINITIONS 


The  technical  terms  defined  herein  supercede  those  given 
in  Volume  I,  AFWAL-TR-86-3017  (2).  New  terms  have  been  added 
and  selected  Volume  I  terms  have  been  revised.  Should  any 
questions  arise,  the  definitions  herein  should  be  used. 
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DEFINITIONS 


1 .  Combined  Least  Square  Sums  Approach  (CLSSA)  -  the 
least  square  subs  for  Individual  fractographic  data  sets  are 
combined  to  estimate  the  EIFSD  parameters  in  a  "global 
sense."  This  approach  is  used  in  conjuncticn  with  the  data 
pooling  philosophy. 

2.  Compatible  Equivalent  Initial  Flaw  Size  Distribution 
Function  -  this  is  a  distribution  function  for  equivalent  in¬ 
itial  flaw  sizes  (EIFS)  which  is  derived  using  a  physically 
meaningful  cumulative  distribution  of  time-to-crack  initia¬ 
tion  (TTCI)  function  and  a  suitable  deterministic  crack 
growth  law. 

3.  Crack  Size  -  is  the  length  of  a  crack  in  a  structur¬ 
al  detail  in  the  direction  of  crack  propagation. 

4.  cumulative  Distribution  of  Service  Time  1F^X> 

-  is  defined  as  the  probability  that  the  service  time  T(x1) 
to  reach  a  crack  size  is  shorter  than'T'.  It  is  equal  to 
the  probability  that  the  crack  size  a (T  )  at  service  life 
T  will  exceed  x^  which  is  simply  the  probability  of  crack 

exceedance ,  i . e . , 

=  ptrCx‘)*T]  * 


5.  Data  Pooling  -  is  a  concept  for  estimating  the  EIFSD 
parameters  using  one  or  more  fractographic  data  seta  in  a 
"global  sense."  A  data  pooling  procedure  is  used  to  increase 
the  sample  size  for  determining  the  EIFSD  parameters. 
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6 .  Deterministic  Crack  Growth  Approach  (DCGA)  -  Crack 
growth  parameters  ara  traatad  aa  deterministic  valuaa  result- 
ing  in  a  aingla  valua  pradiction  for  crack  length. 

7.  Durability  -  is  a  quantitativa  measure  of  the  struc¬ 
tural  resistance  to  fatigue  cracking  under  specified  service 
conditions.  Structural  durability  ia  concerned  with  the  pre¬ 
vention  of  functional  impairments  due  to:  (1)  excessive 
cracking  and  (2)  fuel  leakage/ ligament  breakage.  Excessive 
cracking  is  concerned  with  relatively  small  subcritical  crack 
sizes  (e.g.,  <  0.05")  which  affect  functional  impairment, 
structural  maintenance  requirement  and  life-cycle-coats. 
Such  cracks  may  not  post  an  immediate  safety  problem.  How¬ 
ever,  if  the  structural  details  containing  such  cracks  are 
not  repaired,  economical  repairs  cannot  be  made  when  these 
cracks  exceed  a  limiting  crack  size.  Functional  impairment 
due  to  fuel  leakage/ ligament  breakage  is  typically  concerned 
with  large  through-the-thickness  cracks  (e.g.,  0.50"-0.75H) . 
Although  such  cracks  are  usually  subcritical,  they  affect  the 
residual  strength,  fleet  readiness,  and  may  require  increased 
maintenance  action. 

8.  Durability  Analysis  -  is  concerned  with  quantifying 
the  extent  of  structural  damage  due  to  fatigue  cracking  for 
structural  details  (e.g.,  fastener  hole,  fillet,  cutout,  lug, 
etc.)  as  a  function  of  service  time.  Results  are  used  to  en¬ 
sure  design  compliance  with  Air  Force's  durability  design  re¬ 
quirements  . 

9.  Economic  Life  -  is  that  point  in  time  when  an  air¬ 
craft  structure's  damage  state  due  to  fatigue,  accidental 
damage  and/or  environmental  deterioration  reaches  a  point 
where  operational  readiness  goals  cannot  be  preserved  by 
economically  acceptable  maintenance  action. 
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10.  Economic?  Life  Criteria  -  are  guidelines  and  format* 
for  defining  quantitative  economic  life  requirements  for  air¬ 
craft  structure  to  satisfy  U.  S.  Air  Force  Durability  design 
requirements.  The  economic  life  criterion  provides  the  basis 
for  analytically  and  experimentally  ensuring  design  compli¬ 
ance  of  aircraft  structure  with  durability  design  require¬ 
ments.  Two  recommended  formats  for  economic  life  criteria 


o  probability  of  crack  exceedance 
o  cost  ratio:  repair  cost/ replacement  cost 

11.  Economic  Repair  Limit  -  is  the  maximum  damage  size 
that  can  be  economically  repaired  (e.g.,  repair  0.03"-0.05" 
radial  crack  in  fastener  holes  by  reaming  hole  to  next  size) . 

12.  Equivalent  Initial  Flaw  Size  (EIFS)  -  is  an  artifi¬ 
cial  crack  size  which  results  in  an  actual  crack  size  at  an 
actual  point  in  time  when  the  initial  flaw  is  grown  forward. 
It  is  determined  by  back-extrapolating  fractographic  result*. 
It  has  the  following  characteristics:  (l)  an  EIFS  is  an  ar¬ 
tificial  crack  assumed  to  represent  the  initial  fatigue  qual¬ 
ity  of  a  structural  detail  in  the  as-manufactured  condition 
whatever  the  source  of  fatigue  cracking  may  be,  (2)  no  direct 
relationship  to  actual  initial  flaws  in  fastener  holes  such 
as  scratches,  burrs,  microdefects,  etc.,  and  it  cannot  be 
verified  by  NDI,  (3)  a  universal  crack  shape  in  which  the 
crack  size  is  measured  in  the  direction  of  crack  propagation, 

(4)  it's  in  a  fracture  mechanics  format  but  EIFSs  are  not 
subject  to  linear  elastic  fracture  mechanics  (LEFM)  laws  or 
limitations,  such  as  the  "short  crack  effect"  [e.g.,  41-47], 

(5)  it  depends  on  the  fractographic  data,  the  fractographic 
crack  size  range  for  the  back  extrapolation,  and  the  crack 
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growth  rate  modal  used,  (6}  it  must  be  grown  forward  in  a 
manner  consistent  with  the  basis  for  the  EIFS,  and  (7)  EXFSs 
are  not  unique  -  a  different  set  is  obtained  for  each  crack 
growth  law  used  for  the  back-  extrapolation. 

13.  Equivalent  Initial  Flaw  Size  Distribution  (EIFSD)  - 
is  used  to  represent  the  initial  fatigue  quality  variation  of 
a  structural  detail.  An  EXFS  is  a  random  variable,  and  the 
EIFSO  statistically  describes  the  EIFS  population.  The  EIFSD 
does  not  necessarily  contain  the  "rogue  flaw." 


* 

14.  EXFS  Master  Curve  -  is  a  curve  (e.g.,  equation, 
tabulation  of  a(t)  vs.  t  or  curve  without  prescribed  func¬ 
tional  form)  used  to  determine  the  EXFS  value  at  t»0  corres¬ 
ponding  to  a  given  TTCX  value  at  a  specified  crack  size, 
such  a  curve  is  needed  to  determine  the  EXFS  distribution. 
The  EIFS  master  curve  depends  on  several  factors,  such  as  the 
fractographic  data  base,  the  fractographic  crack  size  range 
used,  the  functional  form  of  the  crack  growth  equation  used 
in  the  curve  fit,  etc.  (Ref.  EIFS). 

15.  Extent  of  Damage  -  is  a  quantitative  measure  cf 
structural  durability  at  a  given  service  time.  For  example, 
the  number  of  structural  details  (e.g.,  fastener  holes,  cut¬ 
outs,  fillets,  etc.)  or  percentage  of  details  exceeding  spe¬ 
cified  crack  size  limits  with  a  certain  probability.  Crack 
length  is  the  fundamental  measure  for  structural  damage.  The 
predicted  extent  of  damage  is  compared  with  the  specified 
economic  life  criterion  for  ensuring  design  compliance  with 
U.  S.  Air  Force  durability  requirements. 


16.  Generic  Sirs  Distribution  -  An  EIFS  distribution  is 
"generic*  if  it  dspsnds  only  on  ths  material  and  manufactur¬ 
ing/fabrication  procsssss.  An  EZFSD  is  not  strictly  "generic" 
bacausa  it  is  basad  on  fractographic  rasulta  which  raflact 
givan  conditions  (e.g.,  load  spactra) .  For  durability  analy¬ 
sis,  an  EIFSD  is  astablishad  using  tha  fractographic  rasults 
for  ona  or  more  data  aats,  and  tha  rasulting  EIFSD  is  justi¬ 
fied  for  a  different  sat  of  conditions. 


17.  Initial  Fatigue  Quality  (IFQ)  -  characterizes  tha 
initial  manufactured  stats  of  a  structural  detail  or  details 
with  respect  to  initial  flaws  in  a  part,  component,  or  air¬ 
frame  prior  to  service.  Actual  initial  flaws  in  a  fastener 
hole  are  typically  random  scratches,  burrs,  microscopic  im¬ 
perfections,  etc.  Such  flaws  are  not  cracks  par  ee  like 
those  associated  with  linear  elastic  fracture  mechanics.  The 
IFQ  is  represented  by  an  equivalent  initial  flaw  size  distri¬ 
bution  (EIFSD) . 

18 .  Probability  of  Crack  Exceedance  (p(i,m  -  refers 
to  the  probability  that  a  crack  in  the  ith  stress  region  will 
exceed  a  specified  crack  size,  at  a  given  service  time, 

T* .  It  can  be  used  to  quantify  the  extent  of  damage  due  to 
fatigue  cracking  in  fastener  holes,  cutouts,  fillets,  lugs, 
etc. 


19*  Reference  Crack  Size  fa0>  -  This  is  the  specified 
crack  size  in  a  detail  used  to  reference  TTCls. 
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20.  Service  Crack  Growth  Master  Curv  (SCGMC)  -  SCGMC 
is  a  curva,  expressed  by  aquation  or  tabulation  of  a(t)  ver- 
sus  t,  usad  to  grow  EIPSs  forward  in  ordar  to  datarmina  tha 
crack  siza  distribution  at  any  sarvioa  time.  Tha  SCGMC  must 
ba  consistant  with  tha  basis  for  tha  EIFS  distribution. 

21.  Sarvica  Time  to  Raach  Any  crack  Siza  x^  -  This  term 
dascribas  tha  time,  T  (x,)  ,  to  raach  any  specified  crack  siza 
x^.  In  this  context ,  tha  crack  siza  x1  can  ba  associatad 
with  aithar  tha  "crack  initiation"  or  tha  "crack  propagation" 
process.  Tha  tima-to-crack-initiation  (TTCI)  term  is  restric¬ 
ted  to  crack  sizes  associatad  with  tha  crack  initiation  pro¬ 
cess,  where  x1  -  af  (reference  crack  siza  for  TTCIs) . 

22.  Statistical  Scaling  -  is  usad  to  account  for  the 
inhomogeneous  fractographic  data,  in  particular  fractographic 
data  associatad  with  tha  largest  flaw  par  specimen  withj. 
holes. 


23.  Stochastic  Crack  Growth  Approach  (SCGA)  -  an  ap¬ 
proach  which  directly  accounts  for  tha  crack  growth  rata  dis¬ 
persion  in  the  durability  analysis. 


24.  structural  Detail  -  is  any  element  in  a  metallic 
structure  susceptible  to  fatigue  cracking  (e.g.,  fastener 
hole,  fillet,  cutout,  lug,  etc.). 

25*  Time-To-crack-Inltiatlon  (TTCI)  -  is  the  time  or 
service  hours  required  to  initiate  a  specified  (observable) 
fatigue  crack  siza,  aQ,  in  a  structural  detail  (with  no  init¬ 
ial  flaws  intentionally  introduced) . 
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26.  TTCI  Lower  Bound  Limit  l  g  )  -  is  a  ainiaua  value 
for  tiae-to-crack  initiation  with  a  reference  crack  size  aQ. 
It  depends  on  the  reference  crck  size  aQ  for  TTCI;  the  larger 
aQ,  the  larger  6. 


27.  Upper  Bound  EIFS  Limit  (xu)  -  defines  the  largest 
EIFS  in  the  initial  fatigue  quality  distribution.  Con¬ 
straints  on  xu  for  fatigue  holes;  largest  EIFS  in  data  set  < 
xu  (e.g. ,  0 . 03"-0 . 05” ) . 
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ACRONYMS 


ADA  «  Advanced  Durability  Analysis 

ASIP  «  Aircraft  Structural  Integrity  Program 


CLSSA 

DADTA 

DCGA 

EIFS  1 

EIFSD 

FHQ 

HEIFS 

IFQ 

LEFM 

LT 

MM 

NDE 

NDI 

NLT 

SCGA 


Combined  Least  Square  Sums  Approach 

Durability  and  Damage  Tolerance  Assessment 

Deterministic  Crack  Growth  Approach 

Equivalent  Initial  Flaw  Size 

Equivalent  Initial  Flaw  Size  Distribution 

Fastener  Hole  Quality 

Homogeneous  EIFS 

Initial  Fatigue  Quality 

Linear  Elastic  Fracture  Mechanics 

Load  Transfer  Through  the  Fastener 

Method  of  Moments 

Non  Destructive  Evaluation 

Non  Destructive  Inspection 

No  Load  Transfer  Through  the  Fastener 

Stochastic  Crack  Growth  Approach 


SCGMC 


Service  crack  growth  master  curve 


SSE 

m 

Sum  Squared  Error 

TSE 

m 

Total  Standard  Error 

TTCI 

m 

Time-to-Craek  initiation 
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list  of  symbols 


a  -  crack  Size 

a0  ■  Reference  crack  siza  for  given  TTCls 

a(0)  -  SZ7S  -  Crack  aiza  at  t*o 

a(t)  -  Crack  siza  at  any  service  time  t 

a(t),  aft^,  a(t2)  -  Crack  siza  at  time  t,  and  t2,  raa- 

pactivaly 

a(T)  "  Crack  siza  at  sarvica  tiaa  T 

a(T)  “  Crack  siza  at  any  sarvica  tiaa  T 

AL,  AU  ■  Lowar  and  uppar  bound  fractographic 

crack  siza,  respectively,  usad  to  da- 
fina  tha  EIFSD  parameters.  Also  usad 
in  conjunction  with  tha  SCGMC  to  da* 
fina  crack  siza  limits  for  tha  snail 
crack  siza  region. 

AU*  -  Uppar  bound  crack  siza  limit  for  tha 

large  crack  size  region 

b,  Q  ■  crack  growth  parameters  in  tha  aquation 

daft)Tg^tj  jb.  Usad  in  conjunction 
dt 

with  tha  IFQ  modal. 

b^,Qj  *  Service  crack  growth  rate  parameters  in 

tha  equation  da/dt  »  Q^a)^  associated 
with  tha  one-segment  DCGA  or  1st  segment 
of  the  two-segment  approach. 


11-1 


tejiQj  •  Service  crack  growth  rate  parameters  in 

the  equation  da/dt  -  Q2(a)b*for  segment 
two  of  the  two-segment  dcga. 

o  •  b  •  l;  used  in  conjunction  with  the  IFQ 

model  when  the  creek  growth  lew, 

de(t)_ jk  is  used  end  b  >  l.o. 
dt 

deft)  ■  Crack  growth  rate  as  a  function  of  time 

dt 


£x(u, 


Probability  density  function  of  X. 


Fft(0) 


(x> 


EIFS  cumulative  distribution  function 
for  a  '‘single  hole  population." 


v> 


(x) 


F*l (0) 
*1 


(Xij) 


F 


e(t) 


(x) 


Fe^(t) (x) 


Cumulative  distribution  of  EIFS  based 
on  the  largest  fatigue  crack  per  test 
specimen  with  l  holes. 

Subscripted  notation  used  for  ra^(o) 

in  conjunction  with  data  pooling,  where: 
j  denotes  the  jth  crack  in  the  ith  data 
set. 


Cumulative  distribution  of  crack  size 
aft)  at  any  service  time  t. 


V. 


Cumulative  distribution  of  crack  size 
a^ft)  at  any  service  time  t  for  the 
largest  fatigue  crack  per  test  specimen 
with  A  holes . 
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TTCI  cumulative  distribution  function 


rT(t) 

Tm  (t) 


F  CT) 

T(»t) 

G(xi;r|x-u) 


X 


Cumulative  distribution  of  minimum 
TTCIs  based  on  the  largest  fatigue 
crack  per  test  specimen  with  JL ^  holes. 


Subscripted  notation  used  for  (t) 

in  conjunction  with  data  pooling,  where: 
j  ■  jth  TTCI  value  in  the  ith  data  set. 


Cumulative  Distribution  of  service  time 
T(Xj  )  to  reach  a  crack  size  x^. 

Initial  flaw  size  corresponding  to  crack  size 
at  time  T  with  X  -  u. 

No.  of  fastener  holes  per  test  specimen. 


L(T),  L(r) 


LT 

m 


M 


J 


“  Total  and  average  number  of  details,  respectively, 
in  the  entire  component  having  a  crack  size  i  Xj^ 
at  any  service  time  TT  • 

■  Load  transfer  through  the  fastener. 

-  Number  of  stress  regions  (or  total 

number  of  fatigue  cracks  in  a  data  set, 
Eqs.  (3-33U3-34)  • 

*  Total  number  of  EIFS  data  sets  used  to 
estimate  the  EIFSD  parameters. 

■  Number  of  TTCI  or  EIFS  values  for  the 
ith  data  set  used  in  conjunction  with 
the  combined  least  square  sums  approach. 
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nu ,n#  ta ,T) 


P(i  ,T) 


Total  and  average  number  of  dataila, 
respectively,  having  a  crack  aisa  ax* 
caading  x^  at  any  service  tlaa  T 

Probability  that  a  datail  in  tha  ith 
atraaa  ragion  will  hava  a  crack  aixa 
>xN  at  tha  aarvica  tine  T 


Crack  growth  rata  parameter  (see  Eq.(  3-6) 
for  the  ith  fractographic  data  set  or 
"pooled  Q"  value.  It  is  used  to  determine 
EIFSs. 


Crack  growth  rata  parameter  (see  Eg. ^ 3-5) 
for  the  jth  fatigue  crack  in  a  fractographic 
data  set. 


tf  tl'  *2 


night  hours  it  t,  t2,  respective¬ 
ly. 


T,  TTCI 


“  Tiaa-to-crack-initiation 


T(xx) 


Service  time  to  reach  any  crack  size  x^. 


A  particular  value  of  X  (lognormal  random 
variable) . 

Crack  size 

Crack  aize  used  for  p(i,T')  predictions  or 
reference  crack  size  for  )(T)  predictions. 
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Upper  bound  limit  for  SZ75 


*u 

X 


Lognormal  random  variable  with  a  median 
of  1.0. 


In  ln(jcu/xij) 


yu<r) 


An  EIFS  value  in  the  EIFSD  corresponding 
to  a  crack  size  x1  at  time  T" in  the  ith 
stress  region. 


Z 


-  Log  X 


n  ) 


Gamma  function 


c,  v 


Empirical  eonatanta  in  the  equations 
^  mctf ,  where  C  *  etreae 


Standard  deviation  of  Z  -  Log  X. 


r 


A  particular  service  time 


Weibull  compatible  shape  and  scale 
EZFSD  parameters,  respectively 


r 
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